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Resumo 



Nesta dissertação, utilizamos métodos da teoria de processos estocásticos para a 
compreensão da turbulência em fluidos. Discutimos o modelo de Kolmogorov para 
turbulência homogénea desenvolvida, resultados analíticos recentes para a equação 
de Burgers. 

Palavras-Chaves:Turbulência; Métodos Estocásticos; Dinâmica de Fluidos 

Áreas do conhecimento: Turbulência, Dinâmica de Fluidos, Sistemas Dinâmicos, 
Processos Estocásticos. 
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Abstract 



In this thesis, we use some methods of theory of stochastic process to approach 

turbulcncc in fluid dynamics. We discuss thc Kolmogorov modcl for fully dcvclopcd 
homogencous turbulence, some recent analytical results to the 1-D Burgers' equation 
approach turbulence. 
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Métodos Matemáticos 
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o objetivo destes dois primeiros capítulos é fornecer o embasamento matemático, 
por mim considerado mínimo para a leitura do corpo principal da tese (capítulos 
3,4 e 5). Embora estes capítulos não sejam imprescindíveis, para uma pessoa com 
experiência em áreas correlatas, ajudam a tornar a tese mais autoconsistente e 
estabelecem os conceitos básicos que norteiam o estudo de turbulência desenvolvida 
na atualidade. 



Capítulo 1 



Sistemas Dinâmicos 



Neste capítulo trataremos de algumas definições e teoremas de sistemas dinâmicos, 
tendo como objetivo final justificar os capítulos posteriores, que encaram o fenómeno 
de turbulência desenvolvida, utilizando ferramentas de processos estocásticos. 

Algumas palavras devem ser ditas sobre o que entendemos como justificar. 
Embora o estudo matemático de sistemas dinâmicos remonte ao começo do século 
com os trabalhos visionários de Poincaré, a aridez deste campo de pesquisa fez com 
que, até esta década, o foco estivesse voltado, quase que exclusivamente, ao estudo 
de sistemas de equações diferenciais ordinárias. Como a passagem dos resultados 
obtidos em EDOs para EDPs permaneça ainda inacabada, este capítulo servirá 
somente para um entendimento qualitativo^ do comportamento das soluções. 

Apesar de o nosso interesse, neste trabalho, estar voltado para sistemas descritos 
por equações diferenciais, existem resultados análogos para mapeamentos. As re- 
ferências para este capítulo são: M, pBl 53, WB, H^, p| . 



1.1 Conceitos Básicos 

Durante quase todo o capítulo estaremos estudando a seguinte equação: 

i = -^x = f{x,fi) (1.1) 

Ou seja, um sistema descrito por uma equação diferencial ( |1.1| ), onde x E U C t, 
/i G 3fí+, com U aberto, í e /i parâmetros, poderíamos pensar que há uma restrição, 

*Como veremos mais adiante a teoria de Sistemas Dinâmicos não é suficiente para o entendi- 
mento da turbulência desenvolvida. 
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pois / não depende explicitamente do tempo, portanto, estaríamos considerando 
somente sistemas autónomos. Todavia, podemos lançar mão do seguinte artifício: 
introduzimos uma variável s, tal que, ^ = 1. Logo, em termos deste novo parâmetro, 
o sistema é escrito como: 

xf=^^fix,s)xl. (1.2) 
Criamos, então, um novo sistema y dado por: 

y = (x, s) 

9{y) = {f{x,s),i), 

obtendo o seguinte sistema : 

= ye3?"x$ft. (1.3) 

Assim, um sistema não autónomo pode ser transformado em autónomo aumentando 
sua dimensão. 

Denotando uma solução com condição inicial Xo{x) no instante Íq por 

x{t) ^ x{t,to,Xo), (1.4) 

podemos ver que esta solução define um fluxo 

(f)t : Ux^+^W" 
(a, í) I— > x{t,Q,a), 

que satisfaz a seguinte equação: 

j^{<P{x,t))U=r^f{<f>{x,T)). (1.5) 

Neste caso, {4>t}'^o forma um semigrupo parametrizado pelo tempo, pois o fluxo 
goza das seguintes propriedades: 

00 = id, 

definindo portanto, uma família de difeomoríismos (f)t :U ^ 
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De um ponto de vista geométrico, x(t, to, xq) é uma curva sobre 3?", parametrizada 
pelo tempo t, tal que, no instante to está em xq G 3?" e sua tangente satisfaz a 



equação diferencial em todos os tempos. Assim, uma solução de equilíbrio 

pode ser vista como um ponto fixo do fiuxo definido por (pt que, por ser a mais 
simples possível, é a primeira a ser estudada. 

Definição 1.1 Chamamos x{t) de ponto fixo ou solução de equilíbrio se f (x) = 0. 

Para entender esta definição, basta substituir em ( |1.1| ) e verificar que a derivada 
temporal desta solução é nula, logo o campo não variará no tempo. Como a perda de 
estabilidade das soluções(conceito que será melhor explicado adiante) é o fenómeno 
que nos leva ao caos, é importante termos definições precisas. Assim, numa primeira 
abordagem, poderíamos exigir que soluções próximas permanecessem próximas para 
qualquer tempo, ou seja: 

Definição 1.2 Uma solução x{t) é dita L-estável (Liapunov estável) se, dado e > 
O, existe 6 > O, tal que para qualquer solução x (t), satisfazendo |x (to) — x (to)| < ô, 
implica \x (t) — X (t)| < e, para todo t > to- 

Definição 1.3 Uma solução x{t) é dita assintoticamente estável, se é L-estável 
e existe a constante real positiva, tal que para toda solução x satisfazendo |ã;(to) — a;(to)| < 
a, temos limí^oo |ã;(t) — x(t)| = 0. 

O que induz as seguintes definições: 

Definição 1.4 Chamamos os pontos fixos assintoticamente estáveis de sorvedouros 
e os L-estáveis, que não são sorvedouros, de centro. 

Uma aproximação um pouco mais física sobre a estabilidade é pensarmos no conceito 
de perturbação. 

Definição 1.5 Seja F G C"^(3fí"), r,k e Z+, k < r, e e > 0. Então, G e é uma 
perturbação de tamanho e, se existe um conjunto compacto K G 3?", tal que, 
F = G Vx G {3?" — K}, e para todo {ii, . . . ,in), com ii + ■ ■ ■ + in = i < k, tivermos 
| L., ^' — iF-G)\\ <e. 
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Definição 1.6 Duas aplicações , F e G são equivalentes ou conju- 
gadas {k < r), se existe um homeomorfismo h G C^, tal que, h o F = G o h. Uma 
equivalência é chamada de equivalência topológica. 

Seja 0í(x, t) um fluxo cujas perturbações e são t). Então, se elas são topologi- 

camente equivalentes, existe um homeomorfismo h, tal que, 4>tl{x, t)oh = hoípl^. (x, í), 
ou seja, existe uma equivalência entre as órbitas geradas por estes dois fluxo. To- 
davia, h não tem que necessariamente preservar a parametrização temporal das 
órbitas, mas quando isto acontece, chamamos o homeomorfismo de conjugação. 

Definição 1.7 Uma aplicação F G (7^(3?) e estruturalmente estável, se existe 
um e > O, tal que, para todas as perturbações de tamanho e de F são topologicamente 
equivalentes. 

A importância destes conceitos reside no fato de que por não serem lineares 
muitos dos sistemas de interesse, necessitamos saber, quando fazemos linearizações, 
quais informações do sistema não linear ainda estão no modelo linearizado, ou seja, 
quais são as propriedades locais. A título de ilustração, consideremos o seguinte 
exemplo: tomemos uma solução arbitrária x, de um sistema com ponto fixo e a 
representemos por: 



onde y é um campo a ser determinado. Agora, substituindo esta expressão na 
equação diferencial e fazendo expansão em série de potência centrada no ponto fixo, 
vem 



X = x{t) + y, 




X = i;{t) +y = í\x{t)) + Vf{x{t))y + 0{\y\^), 



(1.7) 



onde V é operador diferencial no espaço em questão. Mas por (def. |l.lD vem: 





Desprezando termos de segunda ordem em y. Ficamos com 



y = Vf(x{t))y, 




cuja solução é 




(1.10) 



Para analizarmos o sistema linearizado fazemos a seguinte definição: 



Capítulo 1. Sistemas Dinâmicos 



5 



Definição 1.8 Sendo x — x um ponto fixo de x — f{x), x e 3?", x será chamado 
de ponto fixo hiperbólico, se nenhum dos autovalores de T>f{x) tiver parte real 
zero. 

Teorema 1.1 (Hartman-Grobman) Se x é ponto fixo hiperbólico , existe um 
homeomorfismo h definido em alguma vizinhança U{x} C 3?", tal que h{(f){x,t)) = 
e^(^) . Note-se que este homeomorfismo preserva o sentido de órbita e também pode 
ser escolhido de forma que a parametrização temporal seja preservada. 

Teorema 1.2 Se todos os autovalores de T>f{x{t)) têm parte real negativa para um 
tempo suficientemente grande, a solução de equilíbrio de x — x (t) do campo vetorial 
linear é assintoticamente estável. 

Assim sendo, poderíamos nos perguntar se esta condição não é excessivamente 
forte, ou seja, se deixássemos alguns auto- valores serem nulos, será que não teríamos, 
ainda assim, uma solução estável? Infelizmente o método da linearização não nos dá 
a resposta, pois neste caso o comportamento da solução é controlado pelos termos de 
ordem superior em y, na região gerada por esta variedade local central. Mais adiante, 
no entanto, voltaremos a este problema. Para evitar as mazelas da linearização, 
podemos utilizar o seguinte teorema: 

Teorema 1.3 (Estabilidade de Liapunov) Se x é um ponto fixo e y : C/ — > 3? 
uma função definida em alguma vizinhança U de x, tal que: 

1. Vi^) = O e V{x) >Q se, x^ x; 

2. V'{x) < O em U — {x}, então x é estável. 
Além disso, se 

3. V{x) < O em U — {x}, então V é chamada de função de Liapunov. Se 

pudermos tomar U como 3?" , então x é globalmente assintoticamente estável. 

Ao estudar sistemas mecânicos, a energia é sempre uma boa candidata para a 
função de Liapunov, já que é sempre positiva definida e, nos sistemas dissipativos 
sua derivada temporal é negativa. Logo, para sistemas sem termos dissipativos, nos 
quais V{x) > O, este método sempre será inconclusivo. 
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Definição 1.9 Seja -S" C 3?" um conjunto. Então, S é dito conjunto invariante 

sob o campo vetorial x — f{x), se para qualquer Xq & S existir x{t,to,Xo) e S para 
todo í e 3?. Se a assertiva for válida, somente para tempos positivos S será chamado 
de conjunto positivamente invairiante. 

Definição 1.10 Um conjunto invariante S C e chamado de variedade invari- 
ante C^{r > 1), se S tiver a estrutura de uma variedade diferenciável C". 

Com base nestas definições, voltemos a olhar o sistema linearizado. 

Definição 1.11 Seja x um ponto fixo, e o seu sistema linear associado y — Ay, y e 
3?^*; se representamos 3?" como uma soma direta de três subespaços denotados por 

E^, E"^ e E", tais que 

E^ = span{ei, ...Cs}, (l-H) 
= span{es+i,...es+y}, (1-12) 
E" = span{es+i^+i,...,es+i^+c}, (1-13) 

sendo {ei, ...6^}, {cg+i, ...es+u}) {^s+u+i, ■■■■,^s+u+c} os conjuntos de autovetores de 
A{x) com parte real negativa, com parte real positiva e com parte real nula, respec- 
tivamente. Além disso se, 

0í(E") = E", a = s,u,c (1.14) 

para todos os tempos. Então chamaremos E'^ , E'^ e E'^ de subespaços estáveis, 
instáveis e centrais, respectivamente. 

Observe que E^, E'^ e E'^ são subespaços invariantes do fluxo linearizado. Mas, se 
considerarmos o fluxo não linear, eles serão invariantes somente em x. Isto devido 
a contribuição dos termos de ordem maior em y, que promovem uma distorção nos 
autovalores de A conforme nos distanciamos de x. Por exemplo, seja a solução do 
sistema autónomo x — y -\-x., x ponto fixo. Logo, podemos reescrever a EDO do 
sistema como 

y = Vf{x)y + R{y), y^^^, (1.15) 

sendo Riy) = 0(|í/p). Como existe uma transformação linear T, que leva a forma 
linearizada desta equação em uma forma bloco-diagonal 

A, O O W ^ 

O A„ O V , (1.16) 

O Q Ar W W 



V 

w 



( 



\ 



I 
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onde T~^y = {u, v, w) G 3?* ® ® s + u + c = n, Agé uma matriz s x s tendo 
autovalores com parte real negativa, Au é uma matriz nxn com autovalores de parte 
real positiva e é uma matriz c x c com autovalores com parte real nula, aplicando 
esta transformação linear para as coordenadas do problema não linearizado, teremos 

ú = AgU + Rs{u,v,w) 

V = AuV + Ru{u, V , w) (1-17) 
w = AcW + Rc{u,v,w), 

logo, devemos estudar não os subespaços invariantes sobre o ponto fixo, mas va- 
riedades, pois assim, em qualquer ordem de aproximação teremos as propriedades 
desejadas. Estes fatos nos conduzem ao seguinte teorema: 

Teorema 1.4 Suponhamos C"" com r > 2. Então, o ponto fixo {u,v,w) = 

O do sistema possui uma variedade estável local Wi^^^i^O) G , de dimensão s; 
uma variedade instável local Wl^^{0) C , de dimensão u; uma variedade cen- 
tral localWi^J^tí) C C^, de dimensão c. Há uma interseção entre estas variedades 
em {u,v,w) = 0. Estas variedades são tangentes ao campo vetorial linearizado na 
origem(E",a = s,u,c) e, portanto, podem ser por eles representadas localmente. 
Em particular, 

{K V, w)\v = h:{u),w = hUu),VK{0) = VhUO) = 0} 
{{u,v,w)\u = K{v),w = hl{v),VK{0) = VhliO) = 0} (1.18) 
{{u,v,w)\v = K{w),u = hl{w),Vh:{0) = VhliO) = 0}, 

com |f I e \w\ suficientemente pequenos nas expressões para ^^^^^(O), Wl^^{0) e 
^íoc(O)) respectivamente. Vale lembrar que h^iu), hl^{u) , h^{v) , h'^{v) , hl{w) , h'^{w) 
são funções C". Além disso, W(l^{0) e W^/oc(0) ^êm propriedades assintóticas de E^ 
e i?", respectivamente. 

Passemos ao estudo de um outro tipo de solução, chamada de solução periódica. 

Definição 1.12 Uma solução de um sistema autónomo através do ponto xq será 
chamada de periódica de período T, se existir T = mma{a;x{t,xo) = x{t + 
a, Xq), a > 0}, tal que x{t, xq) = x{t + T, xq) para todo t G 3fí. 



wtoM = 
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Considerando um campo vetorial bidimensional sujeito à equação abaixo 

X = f{x,y) (1.19) 

y = 9{x,y), 

com {x,y) G 3Í^, então é verdade que : 

Teorema 1.5 (Critério de Bendixson) Se numa região simplesmente conexa D C 
3fí^ a expressão ff + || não for identicamente zero e não mudar de sinal, então 
não possui órbita fechada, inteiramente contida em D. 



Teorema 1.6 Seja B[x, y) contínua com derivada continua em uma região simples- 
mente conexa não for identicamente zero e não mudar de 
sinal em D, então não possui órbita fechada inteiramente contida em D. 



Existe um conjunto de resultados chamados de Teoria de índices para a análise 
de fluxos bidimensionais, que pode ser resumida da seguinte forma: 

Seja r qualquer órbita fechada num plano que não contenha nenhum ponto fixo 
do campo. Movendo-se, por exemplo, no sentido horário (positivo) por F, os vetores 
que representam o campo rotacionam-se, sendo essa rotação de 27ck para uma volta 
completa em T, e k um número inteiro. Este k é chamado de índice de F. O índice 
de um circuito fechado, não contendo ponto fixo, pode ser calculado pela integração 
do ângulo dos vetores em cada ponto do circuito F, medido em relação a algum 
sistema de coordenadas escolhido. Assim, o valor de é dado por 

E a classificação do ponto em questão é dada pelo seguinte teorema: 

Teorema 1.7 (Teorema dos índices) As seguintes propriedades são válidas: 

1. O índice de um sorvedouro, uma fonte ou um centro é +1; 

2. O índice de um ponto de sela hiperbólico é -1; 

3. O índice de uma órbita fechada é +1; 

4- O índice de uma curva fechada, não contendo qualquer ponto fixo, é 0; 
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5. O índice de uma curva fechada é igual à soma dos índices dos pontos fixos no 
seu interior. 

Corolário 1.8 Dentro de qualquer órbita fechada 7 deve haver, no mínimo, um 
ponto fixo. Se existe um, então deve haver uma fonte, um sorvedouro ou um centro. 
Se todos os pontos fixos em 7 são hiperbólicos, então deve haver um número ímpar 
2n + 1, dos quais n + 1 são fontes, sorvedouros ou centros. 

Para discutir o comportamento assintótico das soluções, cabem as seguintes 
definições: 

Definição 1.13 Um ponto e 3?" é chamado a;-limite de x E W^, se existe uma 
seqúência {ti} ,ti t 00, tal que, lim (x, íj) = Xi_^, sendo o conjunto de todos os 
pontos cu -limite denotado por cu (x). 

Definição 1.14 Um ponto Xq, £ 3?" e chamado a-limite de x & se existe uma 
seqúência {íj},íj t —00, tal que, lim0(a;,íi) = Xa, sendo o conjunto de todos os 
pontos a-limite denotado por a (x). 

Definição 1.15 Um ponto xq é dito não-errante se, para todo aberto U, tal que 
Xq & U existe pelo menos um t ^ O tal que (0 {U, t) ílU) ^ 0. 

Da definição anterior, segue que pontos fixos e órbitas periódicas são pontos não- 
err antes. 

Um conceito que se mostrou útil nos estudos de Sistemas Dinânicos, ao contrário 
dos conceitos de cu, a-limite, foi a idéia de conjunto hiperbólico que conjuntamente 
com o conceito de conjunto não-errante formam a base para definir atratores estra- 
nhos e portanto para formar um cenário de transição à turbulência. 

Definição 1.16 Dizemos que um conjunto compacto e invariante A é hiperbóUco 

para o fluxo (pt : M ^ M, se existe C>0;0<A<1; e seu fibrado tangente T\M 
pode ser continuamente decomposto como a seguinte soma direta: T\M — ® E^, 
onde, para todo x e A e todo t, temos: 

yv G E^,^ II (f).tv II < CA*|| V II; 
Wv e E^,^ II (ptv II < CA*|| V \\. 
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Ou seja, o espaço tangente de cada ponto de um conjunto hiperbólico pode ser 
decomposto como a soma direta de um subespaço estável com um instável. 

Definição 1.17 Um conjunto invariante fechado A C 3?" e chamado conjunto 
atrator, se existe alguma vizinhança U de A, tal que 



Segue da definição que o conjunto A — r\t>o<Pit, M) é um conjunto atrator. 

Definição 1.18 Se B = Ut<o<f) {U,t) , onde U é tal que \ÍTat^oo4'{U,t) = A, onde 
A é o conjunto atrator. Então, B é chamado de bacia de atração de A. 

Definição 1.19 Um conjunto fechado e conexo M é chamado de região de apri- 
sionamento, se (f){t, M) G M,yt> 0. 

Definição 1.20 Um conjunto invariante fechado A é topologicamente transi- 
tivo se, para quaisquer dois abertos U,V (Z A, 3í e 3?, tal que ^{t, C/) fl F ^ 0. 

Definição 1.21 Um conjunto A é chamado de atrator, se A é um conjunto atrator, 
e além disto, for topologicamente transitivo. 

Definição 1.22 Chamamos um fiuxo 4>t{x,t) de fluxo tipo Axioma A se o con- 
junto dos pontos não-errantes íl de (pt goza das seguintes propriedades: compacto, 
hiperbólico e os pontos fixos e as órbitas periódicas são densos em fl. 

Definição 1.23 Chamamos um conjunto U de tipo Axioma A se goza das seguintes 
propriedades: invariante, hiperbólico e os pontos fixos e as órbitas periódicas são 
densos em U e (f) é topologicamente transitivo sobre U. 

Definição 1.24 Chamamos um conjunto A tipo Axioma A de atrator tipo Axi- 
oma A, se existe U, tal que A = r\t>Q(f)t{U) . 

Proposição 1.1 Seja 4> uitT' fluxo e M um conjunto compacto invariante por (j). 
Então, para todo p & M vale: 






1. u{p) ^ 0; 
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2. u!{p) é fechado; 

3. ui{p) é invariante sobre o fluxo; 

4- uj{p) é uma união de órbitas fechadas; 
5. cuip) é convexo. 
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Lema 1.1 Se cu (p) não tem pontos fixos, então cu (p) é uma órbita fechada. 

Lema 1.2 Sejam os pontos fixos Pi ^ P2 ^ ^(p); dentro de uma região de apri- 
sionamento. Então existe, no máximo, uma órbita j G cuíj)), tal que a (7) = Pi e 

^(7) =P2- 

Para fluxos bidimensionais, podemos classificar completamente toda a família de 
conjuntos cu{x) associada aos pontos de um conjunto invariante, pelo seguinte teo- 
rema: 

Teorema 1.9 (Poincaré-Bendixson) Seja p & M d , sendo M uma região 
positivamente invariante para o campo vetorial contendo um número finito de pontos 
fixos. Então, uma das afirmações abaixo, mutuamente excludentes, é verdadeira: 

1. cu (p) é um ponto fixo; 

2. cu (p) é uma órbita fechada; 

3. cu (p) consiste em um número infinito de pontos fixos {Pi}°li e órbitas 7 com 
«(7) =Pi eí^(7) =Pj- 

Ou seja, os atratores de um fluxo bidimensional só podem ser de três tipos: 

(a) pontos fixos; 

(b ) órbitas fechadas; 

(c) união de ponto fixos e órbitas fechadas, geradas pelas intersecções das 
variedades estáveis e instáveis destes pontos fixos, ou seja, órbitas de 
pontos fixos hiperbólicos. 

Teorema 1.10 (Hopf) Se o sistema x = f^ (x), x e W^, e 3? possui um ponto 
fixo (xo,/Uo), tal que as seguintes propriedades são verificadas: 

Hl Dxf/XQ {xq) tem um par de autovalores imaginários, e nenhum outro autovalor 
com parte real nula. 
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Então (Hl), implica que existe uma curva suave de pontos fixos comXfj,^ = xq, 



então, existe, uma única variedade central tridimensional passando por (xo,/io) G 
3?" X 3? e um sistema de coordenadas suaves, tal que a expansão de Taylor até terceira 
ordem sobre a variedade central é dada por 



Portanto temos três casos: 

Se a O, existe uma superfície de soluções periódicas na variedade central que pos- 
sui tangenciamento quadrático com os autoespaços de A^^,, e A* , dado (com 



Se a < O, então estas soluções periódicas são ciclos limites estáveis. 
Se a > O, então estas soluções são repulsivas. 

Teorema 1.11 (Peixoto) Um campo vetorial C' sobre uma variedade bidi- 
mensional diferenciável compacta é estruturalmente estável, se e somente se, 

1. O número de pontos fixos e órbitas periódicas é finito, sendo todos hiperbólicos; 

2. Não existe órbita conectando pontos de sela; 

3. O conjunto não-errante é formado por pontos fixos e órbitas periódicas. 

Além do mais, se é orientável, o conjunto dos campos estruturalmente estáveis 
é aberto e denso em iy'"(M^)[|. 

Definição 1.25 Seja J^t (x) uma superfície transversal a x {t,xo), órbita periódica 
de período T, em um ponto x = x{ti,xo) para algum ti. Chamamos (x) de 
secção de Poincaré. 

Definição 1.26 Chamamos de apHcação de Poincaré, a aplicação P : V d 
Ht {x) Tit{x), tal que para todo x E V, P (x) = (j){x,T), onde T^t {x) é a 
secção de Poincaré do fluxo (p, e V C T^t (x) um subconjunto aberto. 

espaço de Hilbert dos campos C sobre Af^. 




X = 



y = 




correção até segunda ordem) pelo parabolóide fi 
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Para ilustrar estas definições, vamos a um exemplo simples. Seja um sistema 
periódico no tempo, i.e., /(x, t) = /(x, í + 27r), podendo ser facilmente reduzido 

a 

X = f{x,9) (1.22) 
ê = u {x,e)eW'xS\ (1.23) 

Portanto, seu fluxo será (f){x,t) = {x{t),d{t)), com 6{t) = uj{t) + 6Q, {mod27T). Defin- 
imos uma seção transversal S,- {9q) = {{x,9) G 3?" x 5*^; 9 = 9q E [0,27r]}, que é a 
seção de Poincaré, sendo que neste exemplo, em particular, ela é global, exceto o 
caso a; = 0. Da definição de aplicação de Poincaré segue 

ou seja, a aplicação de Poincaré é uma geratriz arbitrária do cilindro infinito S^. 



Definição 1.27 Seja A um espaço métrico compacto, N{r,A) a menor cobertura 
de A com abertos de raio r. Chamamos de dimensão de capacidade, e denotamos 
por dimx A, a quantidade 

iogN(r,A) 

dimií A = — lim sup . (1.25) 

r^o log r 



Definição 1.28 Seja A conjunto não-vazio, a uma cobertura finita, cujo k-ésimo 
conjunto cr^ tem diâmetro dk = dimcTfc < r. Dado a positivo, e (A) = info- {dk)" , 
definimos a dimensão Hausdorff de A, como 

dimjíA = sup|a;limm°(A) > ol. (1.26) 



Teorema 1.12 (Sinai-Ruelle-Bowen) Se 0í : 3?" — > é um fluxo classe e 
possui um atrator hiperbólico A. Então, existe uma única medida n, com suporte 
em A. Então, dada uma g contínua q.t.p.fX]^ vale: 

lim ^ r g{(t)t{x))dt= í gdf^ . (1.27) 

T^oo I JtQ Ja 
í Abreviação para quase por toda parte em relação à medida de Lebesgue. 
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1.2 Um Mecanismo Matemático para a Transição à Tur- 
bulência 

O primeiro cenário para a transição à turbulência que existiu, foi o proposto Landau 
|4|. Neste cenário, o fluxo vai passando por bifurcações de Hopf sucessivas, com o 



aumento do Número de Reynolds. Ou seja, o fluxo migraria de um n-tóro para um 
(n+l)-tóro, numa série infinita! Embora este cenário concorde com a intuição, ele 
não é correto, pois experimentalmente e em simulações numéricas [0, 0, após o 
3-tóro, o fluxo já se torna turbulento, e não somente aperiódicofl[^ . 



Tentemos agora fazer uma reconciliação dos conceitos expostos acima com a 
situação encontrada na equação de Navier-Stokes bidimensional, já que para a tridi- 



mensional nem mesmo temos provas de existência ||52|| . Identificaremos elementos da 
Teoria de Sistemas Dinâmicos na equação de Navier-Stokes e ordenaremos os resul- 
tados acima, mostrando como estes constroem quadros de transição à Turbulência. 
Escrevendo a equação de Navier-Stokes como 

^x = Xfí,(x), (1.28) 

para Re = O, o sistema tem uma solução x = 0. Consideremos agora, o jacobiano 
desta solução 



4 = ^ (0) (1.29) 



que tem parte real negativa, ou seja, x = O é um ponto fixo atrativo. Assim, o 
determinante do jacobiano é não nulo. Então, pelo Teorema da Função Implícita, 
existe {t), chamada solução estacionária, dependendo continuamente de Re, tal 
que 

XRei^Re)=0, (1.30) 

OU seja, para todo Re existirá uma solução de ponto fixo(solução estacionária). 
Mas, por continuidade, o jacobiano de ^Re, denotado por Jl [^Re], para um número 
de Reynolds suficientemente pequeno, deve ter parte real negativa. 



§Ondc consideramos um número de Reynolds suficientemente grande(10^ 
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Ao aumentarmos o Número de Reynolds, a parte real dos autovalores pode se 
anular, fazendo com que os pontos fixos se unam ou desapareçam. Ou que suces- 
sivos pares de autovalores complexos possam cruzar o eixo imaginário. Logo pelo 
teorema de Hopf o que era um ponto fixo(solução estacionária), passa a ser uma 
órbita periódica(solução periódica), sendo sua amplitude proporcional ao Número 
de Reynolds. Todavia, à medida que aumentamos Re, podemos ter uma outra bi- 
furcação de Hopf ou transformação de uma órbita periódica em uma quase-periódica. 
Sendo assim, pelo Teorema de Peixoto, a vizinhança desta órbita está em um con- 
junto que é o complemento de um subconjunto denso e aberto do espaço de Hilbert 
dos campos C" sobre o tóro. Este conjunto, em diversas situações, é um atrator, 
sendo este conjunto que defi ne o chamado cenário de transição à Turbulência. 

O mais conhecido cenário, por ser o mais antigo(1971,1978), é o chamado Cenário 
de Ruelle-Takens, embasado no seguinte teorema: 

Teorema 1.13 (Newhouse- Ruelle-Takens) Seja v um campo vetorial constante 
sobre um tóro T" = SÇ^/Z". Se n > 3, toda vizinhança de v contém um campo 
vetorial v com um atrator estranho tipo Axioma A. Se n > A, podemos tomar C°° 
em vez de C"^ . 

Como o teorema exige n > 3, então o sistema deve passar primeiro por três 
bifurcações de Hopf antes de encontrar o atrator estranho[] tipo Axioma A. 

O cenário de Feigenbaum consiste em um processo de duplicação infinito de 
períodos da órbita em um intervalo finito de variação do Número de Reynolds, 
sendo este atrator aperiódico, ou seja, composto por órbitas estáveis de período 
2°°. Assim, o sistema migra para a fase turbulenta por intermédio de bifurcações, 
como as do tipo pitchfork até que o Número de Reynolds atinja seu valor crítico e 
o sistema se estabeleça em um atrator aperiódico e ergódico. 

Como fica patente pela exposição acima, a aproximação, através da Teoria de 
Sistemas Dinâmicos, não nos traz informações sobre a estrutura espacial do fluxo. 
Por exemplo, não nos diz quais são os eventos físicos associados à mudança das 

propriedades do espaço de Hilbert das soluções, existência das regiões inerciais, leis 

^ Um atrator tipo Axioma A é chamado de estranho, se apresenta sensibiUdade em relação as 
condições iniciais, i.e., a distância de duas órbitas cresce exponencialmente no tempo. 
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de potência. Logo, se desejarmos caracterizar a estrutura espacial das soluções 
assintóticas, é muito útil lançar mão sobre a hipótese de ergodicidade, via Teorema 
Sinai-Ruelle-Bowen. Partindo desta hipótese, que estes atratores são ergódigos, 
podemos trocar os cálculos espaço-temporais por médias sobre o atrator e vice- 
versa, ou seja, estaremos fazendo uma análise estatística da Turbulência (ver Ruelle 
emg). 



Capítulo 2 



Processos Estocásticos 



Este capítulo vai no mesmo espírito do anterior, assim que, para as pessoas sem 



contato prévio com o assunto, aconselhamos as seguintes referências: [0, 0, |3§, 



2.1 Conceitos Básicos 

Definição 2.1 Uma tripla {VL,T,P) é chamada de espaço de probabilidade, se 

VL é um conjunto não-vazio, T uma o-álgehra minimal sobre os borelianos de Q, e 
P : Q ^ E uma medida de probabilidade sobre o espaço linear E. 



Definição 2.2 Uma variável aleatória x sobre o espaço de probabilidade{Q, JF, P) 
é uma função x : Q ^ E mensurável em relação a T . 



Definição 2.3 Um Processo Estocástico é uma coleção de variáveis aleatórias 
X definidas em um mesmo espaço de probabilidade (Q,J-',P), indexados por ele- 
mentos de um conjunto de parâmetros T. Quando T = Z, , também chamamos o 
processo de sequência aleatória, se T = Si" chamamos de campo aleatório ou 
função aleatória. 



Definição 2.4 Se X é um processo estocástico, então para t G T, tu G íl,Xt{uj) = 
x{t,-) : Q —>■ E, é mensurável em relação a T; por outro lado Xi^{t) = x{-,uj) : 
T ^ E é contínua em relação à T , sendo chamada de trajetória ou realização 
do processo. 
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Definição 2.5 Seja {ti}^^i G T, uma sequência nos índices e C G B{E"') um bore- 
liano no espaço linear E'^. Definimos a distribuição finita por Pt^ , t^{C) = P{{uj G 

Q]{x{ti,Uj),...,x{tn,Lj))eC}) . 

Assim, a distribuição finita é uma medida de probabilidade sobre B{E^), já que 
x{ti, ■) é jF-mensurável para cada ti G T. Há um teorema de Kolmogorov que 
garante a recíproca da afirmação anterior. 

Teorema 2.1 (Kolmogorov) Dado um conjunto T e uma família de medidas {Pti...t„} 
sobre um espaço linear E satisfazendo: 



1- Pu...t,...tÁC) = Pt,t,...t,...tÁC), vara todo C G B{E^),l < j < n, ou seja, 



invariância por permutação em T. 
2. Pt,...t^...t„{yi,...yj,oo,...,Qo) = Pt,...t^{yi, ...yj). 

Então, existe, pelo menos um processo aleatório X , definido sobre um espaço {íl, JF, P), 
tal que as medidas {Pti...tn} ^ão funções de distribuição finita do processo X . 

Portanto, não precisamos nos preocupar, em primeira análise, com a medida no 
espaço de probabilidade, pois as distribuições induzidas no espaço linear, no nosso 
caso o 3?", já nos dão informações suficientes sobre o processo estocástico em questão. 

Definição 2.6 Definimos o n-ésimo momento da variável aleatória x por: 



pois, Px{k) é a distribuição finita, não a densidade de probabilidade. Vale lem- 
brar aqui que a derivada tomada em (|2.2|) , é a derivada no sentido do teorema 
de Radon-Nikodym, ou seja, ^P^^k) é uma distribuição, e não é necessariamente 
absolutamente contínua em relação à medida de Lebesgue. 




(2.1) 



Logo, para E = ^ teremos: 




(2.2) 



Definição 2.7 Se um processo é tal que para todo elemento de T 
então o chamamos de processo de segunda ordem. 
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Observe que podemos construir classes de funções iguais em quase todos os pon- 
tos e quadrado integrávelQ Essas classes, vistas como elementos de um conjunto, 
podem ser consideradas como um espaço de Hilbert L2(Í2,JF, P), quando definimos 
a norma e o produto interno pelas relações: 

(x, y) = E{xy*), \\x\\ = ^/E{\xf). (2.3) 



Definição 2.8 Definimos a covariância, chamada em turbulência de correlação, 

por K{s, t) = E{xsXl) = {xs, x^) com s,t E T. 



Agora que já temos os conceitos básicos, passemos para a próxima etapa, que é a 
construção do cálculo para processos estocásticos. 

A abordagem escolhida é a que nos parece mais intuitiva, por ser mais parecida 
formalmente com o cálculo usual. Adotaremos, doravante, T = 3?, e assumiremos 
que todos os processos estocásticos são de segunda ordem, ou seja, tem covariância 
finita. Com isto esperamos facilitar o entendimento. Mas, cabe ressaltar que existem 



trabalhos, nos quais sao considerados processos mais gerais p|, |4J para a descrição 
dos vórtices do campo de velocidade de um fluxo turbulento. 



2.2 Cálculo Diferencial 



Diremos que liuit^to Xt = y,se y E L2 e hmt^to E{\xt - y\ 
que o limite existe se ele converge na média quadrática. 
Assim definimos 



d 

—Xt = Xft, 

dt 



se xft E L2 e 



lim 



Xt+h - Xt 



- x\ 



o, ou seja, diremos 



(2.4) 



(2.5) 



Obviamente diremos que Xt é diferenciável, se a derivada existe, neste sentido. 



*Funções que diferem sobre um conjunto de medida nula são identificadas e geram uma classe. 
As classes, produzidas sob esta identificação, geram um espaço que goza da propriedade de todos 
os seus elementos (classes de equivalência), serem quadrado-integráveis. 
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Como os processos de segunda ordem sao caracterizados, a menos de uma isome- 
tria no espaço de Hilbert, pela covariância, a proposição abaixo é muito útil. 

Proposição 2.1 Dado um processo estocástico de segunda ordem X valem as seguintes 
propriedades: 

1. SeYmit^toXt existe-^^ lims^t^to K{s,t) existe. 

2. Se {xt} é contínuo em t = to <í==^ K{s,t) é contínuo em {to,to). 

3. xft existe, se somente se, 

lim — {K(t + h,t + k) -K(t + h, t) - K(t, t + k) + K(t, t)} (2.6) 

h,k^O hk 

existe. 

4. {xt} é continuamente diferenciável em (a, b) <í=^ -^^K(s, t) existe e é contínua 
Ws, t G (a, b). 

Definição 2.9 Um processo é dito estacionário se suas distribuições finitas são 
invariantes sob translação f\, ou seja, se para todo n fixo e toda sequência {tj^}^^^ G 
^, ehe^ tivermos Pu...tSC) = Pu+h...t^+h{C),C e B{E^). 

Definição 2.10 Um processo X é dito estacionário num sentido amplo, se a 

função covariância é invariante sob translação, ou seja, K{s,t) = K{s + h,t + h) . 

Definição 2.11 Dizemos que o processo é de incrementos estacionários, se para 
todo n fixo e toda sequência {t/tj^^i G 3fí, e G 3? a distribuição conjunta dos 
incrementos {xt2 — a^tj, {xt„ — a;í„_i) é invariante pela troca t^ ti + h. 

Definição 2.12 Chamamos um processo W de movimento browniano, se for 

de incrementos estacionários, independentes, contínuo na média e de quadrado in- 
tegrável. 

Teorema 2.2 Todo movimento browniano W com = O , centrado na média, 
goza das seguintes propriedades: 

T Normalmente, quando a família de parâmetros T = SR" representa um processo no espaço 
físico, o chamamos de homogéneo; se representa o tempo chamamos de estacionário. 
■I-Os matemáticos chamam- no de processo de Wiener 



Capítulo 2. Processos Estocásticos 



21 



1. E{Wt)= ^^t,^^ = E{Wl). 

2. E{Wf) = aH,a^ = E{Wl). 

3. K{s,t) = a^mm{s,t). 

4- W é continuo, mas não diferenciável(em relação ao tempo). 
Prova 



1. Observe que: E{Wt+s) = E{Wt)+E{Wt+s-Wt) = E{Wt) + E{Ws), como W é 
contínuo na média, a expectativa deve ser contínua, como E{Wo) = O teremos 

!•)• 

2. Por analogia com o item 1.), vem: E{Wl^,) = E{W^) + E [{Wt+s - Wtf] = 
E{Wt^) +E [{W, - Wt^tY] = E{Wt^) +E{W^), pelas mesmas razões de 1.) 
segue o resultado 2.). 

3. Por ser K{s,t) = E{WsWt) = \ {^(W^^) + Ei^^l) - E{\Wt - W^f)] = 
^{t + s-\t-s\}. 



4. Por ( p.6|) , a derivada do processo só existe se o limite equivalente de covariância 
existe. Como é fácil ver por 3.), este limite não existe, logo W não é difer- 
enciável. 



2.3 Cálculo Integral 

Seja um processo estocástico X, assumindo valores em um espaço de Banach V 
arbitrário, P uma partição de [a,b] G 3?, S{P) = max{(íj_|_i — ti);ti+i,ti G 3?}, o 
maior intervalo da partição P, Q G um ponto do i-ésimo intervalo, onde 

0<i<n-l. 

Definição 2.13 Definimos soma de Riemann da função {xt} o vetor: 

n-l 

Sp = J2^(Áti+i-ti). (2.7) 

i=0 
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Definição 2.14 Se o limite da soma de Riemann, quando Ô{P) — > O, existe, dize- 
mos que {xf} é Riemann integrável em [a, b] e escrevemos: 



rb 

lim Sp= Xtdt. (2.8) 

P)^0 Ja 



SiP) 



Proposição 2.2 A Ja Xtdt existe, se e somente se,J^ds dtK{s,t) existe no sentido 
de Riemann e vale a igualdade: 



pb rb rb 

ds dtK{s, í) = / Xtdt 

Ja Ja Ja 



(2.9) 



Proposição 2.3 Seja y quadrado-integrável. Então, o produto interno e a inte- 
gração comutam, ou seja, vale a seguinte igualdade: 



y, f dtxA = í'dt{y,xt). (2.10) 

Ja I Ja 



Fica claro que podemos estender vários resultados clássicos para esta integral, 
tais como o Teorema Fundamental do Cálculo, integrais impróprias e integração por 
partes. Existe outra formulação mais sofisticada, que se assemelha a de Lebesque, 
que nos permite fazer integrais sobre processos estocásticos, embora seja mais ab- 
strata. 

2.3.1 Integração estocástica 

Como já vimos, as variáveis aleatórias induzem uma medida sobre a álgebra gerada 
pelos subconjuntos de seu contradomínio, chamada função de distribuição finita. 
Donde, surge a seguinte pergunta: Quais as condições que um processo estocástico 
deve satisfazer para induzir uma medida no espaço das funções de T — > com T 
conjunto de parâmetros e V um espaço de Banach? 

Buscamos definir Jrpfdxt,\ff £ L2{T,dF) tal que, dF seja uma medida que 
dependa de um processo esocástico. Ou seja, desejamos construir uma integral es- 
tocástica que nos permita saber o valor médio de funções que dependem de processos 
estocásticos, por exemplo: valores esperados e representação espectral de processos 
estocásticos. 



Capítulo 2. Processos Estocásticos 



23 



Definição 2.15 Dizemos que uma função aleatória possui incrementos ortogo- 
nais, se V(a, 6), (c, d) G T; (a, b) fl (c, d) = tivermos {xb — x^, — Xc) = O 



Teorema 2.3 S'e7'a {xt} uma função aleatória de segunda ordem com incrementos 
ortogonais. Então, existe F : T função não decrescente real, tal que, F é uma 

medida sobre a reta induzida pelo processo {xt}. Sendo F definida por: 



F{t) = E{\xt-Xt,\^},t>to 
F{t) = -E{\xt-Xtf},t<to. 



[2.11) 



Por ( p.ll|) fica claro que: Xt é contínuo em t, se somente se, F{t) é contínua em t . 
Agora estamos aptos a construir a integração estocástica. 



n-l 

Definição 2.16 Seja f{t) = J2 (^i^ (ti u+i){t) uma função escada onde cada ti é ponto 
de continuidade de F. Definimos a integral estocástica de f pelo processo {xt} 
por: 

» n—l 

/ fdxt = Ci{xt^^^ - XíJ. (2.12) 

O ponto importante é que a correspondência fdxt ^ f é uma isometria entre 
o espaço das variáveis aleatórias da forma dada pela equação acima e o conjunto 
das funções- escada consideradas como elementos de L2{T,dF). Mostremos esta 
afirmação por intermédio de um cálculo simples: 

n—l n—l 







j^fdxt 









1=0 j=0 



n-l 



tb-í ^ 2^ n—l r. 

EIqI^í? k,,-xj = Y.hf[Fit.+i)~Fit.)]= \fit)\'dF{t). 

i=o ^ ■> i=o 

Como para toda /, Borel mensurável em T, existe uma seqiiência crescente 
de funções- escada^ {fn}, tal que, lim/„ = /, a seqiiência de funções-escada é 
uma seqiiência de Caucliy. Então, usando o Teorema da Completeza para o 
espaço L2{T,dF), temos assegurada a relação / G L2{T,dF). Mas pela isome- 
tria, {Jrp fndxt} é uma seqiiência de Caucliy em ^2(^2, J-", P). Em seguida, aplicando 
o Teorema da Convergência Monótona vem lim Jj, fndxt = Jj- fdxt. Vale a pena 



§0 que está por trás desta argumentação é que o conjunto das funções escadas é denso neste 
espaço. 
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notar que a integração não precisa ser sobre T, mas pode ser sobre qualquer con- 
junto de Borel em T. Desta forma, a integração estocástica gera um novo processo 
estocástico com incrementos ortogonais. 



Parte II 
Turbulência Desenvolvida 
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Capítulo 3 



Teoria de Kolmogorov 



Os modelos para a turbulência são devidos à inexistência de condições que fechem 
o sistema, ou seja, o cálculo da correlação de segunda ordem envolve a de terceira 
ordem, o cálculo da correlação de terceira ordem envolve a de quarta ordem, e 
assim sucessivamente, levando a um sistema indeterminado, infinito, e acoplado de 
equações diferenciais parciais. Um método para atacar esta situação é fazer uma 
hipótese que feche o sistema em alguma ordem. 

O modelo de Kolmogorov [0^ , ^ foi sugerido em 1941 para explicar quais as 



condições suficientes para obtermos a chamada Lei dos 2/3Q e nos dá também 
uma expressão para a correlação de terceira ordem. Buscando um melhor entendi- 
mento, seguiremos por etapas. A obtenção da relação de Kárman-Howarth ||33|| é 



um primeiro passo para o entendimento do problema, pois gera o cenário onde as 
diversas teorias entrecruzam-se e buscam respaldo teórico e experimental 0, De- 
pois explicaremos brevemente as leis experimentais; passaremos para as hipóteses de 



Kolmogorov e suas implicações limítrofes; exporemos os argumentos de Landau |4Õ 
e Kraichnan||18||, contra a universalidade da constante de Kolmogorov. 



*A Lei dos 2/3 diz que a correlação tem um comportamento assintótico tipo r^/^ para um fluxo 
isotrópico e homogéneo 
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3.1 Simetria e Tensores de Correlação 

Consideremos, primeiramente, o caso da correlação de segunda ordem das compo- 
nentes da velocidade entre os pontos P (xi, X2, x^) e P' {x'i,x'2, x'^) definida por: 

E{u,u',,} = E{u ■ ujR^k. (3.1) 

Calculando a correlação entre a componente da velocidade u {P,t), baseada em P 
numa direção arbitrária Oj e a componente da velocidade u' = u {P', t), baseada em 
P' numa direção arbitrária bk, temos: 

E {u ■ u} R{a,b) = E {uiu[} aibk, (3.2) 

mas por ( p.l|) podemos reescrever a relação acima como 

R{a,b) = Rikaibk- (3.3) 

Se impusermos que R{a, b) deva ser invariante sob rotações de corpo rígido[] e 
translações, para executarmos a invariância translacional, basta termos uma de- 
pendência funcional somente na distância relativa, ou seja, em = (x^ — x,), a,, bk- 
Por conseguinte, a pergunta que se impõe é: Qual a função mais geral R (r; a, b) da 
forma de ( p.3|) que é invariante sob rotações arbitrárias dos vetores ri,ai,bk ? 

A princípio poderíamos pensar que qualquer contração dos três vetores(r, a, b) 
pudesse aparecer na função para R{r;a,b). Entretanto, por (^.3]), R{r;a,b) deve 
ser uma função bilinear em a e b, i.e., não devem participar termos como a - a e b-b. 
Portanto R pode ter no máximo a seguinte forma: 

R (r; a, b) = A{r • a)(r ■ 6) + B{a ■ b) + C{rab), (3.4) 

onde A,B,C funções de r^. 

Comparando ( ^.31 ) com a equação acima, é fácil concluir que: 

Rik = A{r'')rirk + ^(r^)^,, + C{r^)e,uri. (3.5) 
tisto quer dizer que as reflexões são deixadas de fora. 
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Agora voltemos nossa atenção para o tensor de correlação de terceira ordem a 
dois pontos, que é definida por: 

E {uiU.u',} = [E{u- u}f^ Ti.kir). (3.6) 

Exigindo homogeneidade (invariância translacional) e invariância rotacional do 
escalar 

T (r; a, b, c) = TijkaibjCk, (3.7) 

onde Oj e bj são vetores baseados em P e baseado em P'. 

As contrações possíveis, já que T (r; a, b, c) tem que ser linear em a, b, c, são: 

r ■ r, r ■ a, r ■ b, r ■ c, a ■ b, a ■ c, b ■ c, (rabc), (abe). (3.8) 

Logo sua forma é: 

T{r]a,b,c) = Ti{r ■ a){r ■ b) {r ■ c) + T2 {r ■ a) {b ■ c) + (3.9) 
+T3 (r ■ 6) (a ■ c) + T4 (r ■ c) (a ■ 6) + T5 {rabc) + Tg (aòc) . 

Por ( ^.6| ) e pela equação acima vem: 

Tyfc(r) = TiTirfk + ^2^^»^^?= + T^^r^Sik + T4rfc(5íj + T^eijkiri + Tgeíj-fc, (3.10) 

onde os Tj/s são função de r^. 

3.2 A Relação de Kármán-Howarth 

Consideremos a equação de Navier-Stokes, tal que, o campo de velocidades e a força 
sejam estatisticamente isotrópicos e homogéneos. 

dtUi + UjdjUi = -dip + fi + ydjdjUi 
diUi = 0. 

Multiplicando a equação acima por Uk{x + r, í) = e tomando a média no espaço 
amostral-f, vem: 

E {u'i^dtUi} + E {u',^UjdjUi} = -E {u'^.dip} + E {u'Ji} + vE {u'^.djdjUi} 

ÍNote que a correlação entre velocidade e pressão é nula, já que é imposta a condição de isotropia 
e incompressibilidade. 
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E {uAu',} + E {u.fli'^d]u'^] = -E {u,d'^p'} + E {u^fl} + vE [u,d'^d'^u'^] . 

Somando as equações acima, e lembrando que para processos estacionários a co- 
variância depende somente da diferença dos argumentos, a saber r, di = —dr- = — 
e os operadores envolvidos comutam com a média, teremos: 
1 1 

-dtE{uiUk} + -dr^ [e {uiu'^u'i^} -E{u[ujUi} 

= E {Uif-} + vdr^dr^E {UiU'k} , 

podendo também ser escrita sob a forma de: 

-dtRik — -drj [Tjki + Tijk] = Sik + udr^Rik, (3-11) 

onde 5"^^ = E{uifí}. 

A equação (|3.11| ) é chamada de Relação de Kármán-Howarth (RKH). Tomando 
r = O na RKH, esta transforma-se em uma equação para a variação de energia 
do fluxo médio. O primeiro termo do lado direito diz que a força injeta energia 
diretamente no fluxo médio; o segundo termo, que a única forma de dissipação é 
através da viscosidade. E interessante chamar a atenção para o fato da RKH não 
ser um sistema fechado, já que para conhecermos a correlação de segunda ordem 
devemos conhecer a de terceira ordem. 

3.2.1 Análise do Caso Isotrópico 

A hipótese de isotropia^] traz uma grande simplificação na obtenção do espectro da 
covariância, que passa a ser equivalente a do espectro unidimensional, e também um 
ganho no entendimento do mecanismo da turbulência. 

Reescrevendo as correlações na forma de ( ^.51 ) e ( |3.10D , e lembrando que no caso 



isotrópico Rik e Tijk têm que ser invariantes sob reflexão, os termos C,Tr,,TQ devem 
ser nulos. Pela incompressibilidade do fluxo vem: 

dr^Rik = dr^Rik = o (3.12) 



§Nesta secção, exigimos que o termo de força fosse escrito como o gradiente de uma função, 
e portanto, encorporado à pressão e posteriormente elimidado da RKH, via homogeneidade e 
isotropia. 



Capítulo 3. Teoria de Kolmogorov 



30 



Reescevendo as funções de Kármán como: 

f-9 



-íl — õ ) -'4 — — , -'2 — -is — 



teremos por ( |3.12| ) e ( |3.13|) que: 



9 
O 



r d 



^{k~h) + -{k-h-2q) 
or r 



+ 



28,. 



1 ^ 



Substituindo em Q e (CT ) vem: 



Rik{r) 
Tijkij) 



Hk 



Definindo as transformadas de Fourier unidimensionais por: 



E{u-u}fir,t) 
F{k,t) 

E{u ■ u}g{r,t) 
G{k,t) 

substituindo em ( p.l4| ) vem: 



lJZodkexp{tkr)F{k,t) 



E{u.u} roo 



JZcdrexp {-ikr) f{r, t) 



è j^oD dk exp{ikr)G{k,t) 



E{U'U} roo 



2 



IZodrexp{-ikr)g{r,t), 



d 



levando-nos à seguinte representação tridimensional da correlação: 



Rij{r, t) 



dk^exp {ik ■ r) ^ij{k,t) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



1 

^ij{k,t) = j^—^ J dr^ exp {-ik ■ r) Rij{r,t), 

onde, obviamente, Rii{r,t) = Rii{—r,t). Definindo E{k,t) = 2nk'^^ii{k,t), podemos 
escrever a soma das componentes diagonais da correlação como: 

sin(A;r) 



Rii{r\t) 



dkE{k,t)- 



kr 



(3.19) 
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igualando (|3.15|) e (p.l9|) obtemos: 



k-^[F + 2G] = -2E. (3.20) 
ok 



Substituindo ( ^.18| ) na expressão acima vem: 



E{k,t) = V3F{k,t) (3.21) 
^ -2['d^-'dk)- 



A equação acima nos diz que, dado o espectro da função de correlação de se- 
gunda ordem unidimensional, podemos obter o espectro tridimensional da correlação 
de segunda ordem. Analogamente, por (|3.15| ), podemos expressar a correlação de 
terceira ordem, por meio da função de correlação unidimensional^] k = ^^^^—^^jr- 
Importante frisar que esta possibilidade existe em virtude da incompressibilidade do 
fluxo, que nos permitiu escrever ( p.l2| ). 

3.2.2 Transformada de Fourier da RKH 

Olhemos para a RKH, desta vez no espaço de Fourier, para entender por intermédio 
de quais mecanismos a turbulência pode ser gerada, mantida e amortecida. A ( |3.11| ) 
para a direção paralela aos pontos em questão é: 

I (ll«ll V) + 2 \\uf' ( P, + 4^) = 2. Il.ll^ (Pj + ~f), (3.22) 

ot ^ ' \or r \ or^ r or / 



lembrando que: / oc (uiu^) e h oc {ulu'i). Substituindo (|3.16| ) vem: 

|— F + Ak'^Hi -sj sHi{s)ds = -2z/PF + A J sF{s)ds (3.23) 

ko ko 

Hi{k) = / h{r)e''''dr. (3.24) 

/CTT J 



Por ( p.21|) obtemos uma equação para o espectro da energia, aplicando o operador 
Vs em (13:231 ). 

d 

—E = -W-2uPE (3.25) 
dt ^ ' 

W = 4k'^V3[Hi{k)]. (3.26) 



^Estas funções também são chamadas de funções de Kármán-Howarth ou correlações de 
Kármán-Howath. 
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A equação (|3.25|) nos diz que a perda de energia dos turbilhões de escala l/k 
é causado pela dissipação direta 2vk'^E e pela transferência de uma escala para 
outra W . Como W é essencialmente a transformada de Fourier da correlação tripla, 
esta é a responsável pelo transporte de energia de uma escala para outra, sendo a 
correlação tripla herança do termo não linear da equação. Este fato corrobora o 
modelo fenomenológico de Prandtl, onde o mecanismo responsável pelo transporte 
de energia é a colisão dos turbilhões. 

Nos modelos fenomenológicos 0, |2^, ^ a função W é modelada com o 



intuito de ajustar os resultados experimentais. Mas, como estes modelos assumem 
independência estatística entre os turbilhões da região energética e da região dis- 
sipativa, o que é uma aproximação grosseira; e não levam em consideração que os 
turbilhões são bem localizados tanto espacialmente como temporalmente; não são 
capazes de prover um entendimento da dinâmica da turbulência, muito embora, 
dêem resultados satisfatórios na região inercial. 



3.3 Leis Experimentais e Escalas em Turbulência 

O objetivo desta secção é identificar escalas de movimento nos fluidos e comentar 
brevemente as leis experimentais da turbulência. 

Consideremos um fluxo com uma dimensão característica L, velocidade U e vis- 
cosidade cinemática v. Geralmente tomamos U como a raiz quadrada da média 
espacial e amostrai e L como sendo l/fco, onde é ponto de máximo global da 
transformada de Fourier da covariância da velocidade a dois pontos. Neste caso, L é 
chamado de escala energético^, pois os turbilhões que contribuem com uma parcela 
maior, geralmente, estão nessa vizinhança. Chamamos de tempo característico de 
advecçãof^ a 

ta = ^, (3.27) 

pois este dá a ordem de grandeza do tempo gasto por um turbilhão para percorrer 

"Esta é uma tradução totalmente livre do termo: energy-containing scale corrente na literatura 
em língua inglesa. 

**Há aqui um abuso dc linguagem pois o termo advecção refere-se ao movimento de um escalar 
passivo sobre fluido, caso este não considerado na definição. 
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uma distância L no fluxo. Definimos também tempo de difusão 

td = — , (3.28) 

que nos dá a escala de tempo dos processos de difusão viscosa no fluido. 

A existência destes tempos característicos nos sugere a interpretação do Número 
de Reynolds como uma razão entre eles. 

Re = -. (3.29) 

ta 

Como Re é o parâmetro que controla a turbulência, fica patente que a con- 
corrência entre a difusão, advecção e força é o que gera a turbulência, pois a difusão 
tende a diminuir o valor dos expoentes de Liapunov do fluxo, enquanto que o termo 
advectivo acopla os modos normais e a força compensa a perda, fazendo com que 
o sistema não tenda para um ponto fixo. O que qualitativamente concorda com a 
análise feita para EDO com órbitas homoclínicas|^, onde a intersecção entre as 



variedades instáveis e estáveis, em uma configuração tipo ferradura de Smale, é o 
mecanismo gerador de caos. 

De um ponto de vista mais físico, esperaríamos que os fenómenos difusivos fossem 
de pequena importância para a análise dinâmica da turbulência. Todavia, o termo 
viscoso é responsável pela dissipação de energia cinética do fluido em forma de 
calor, fenómeno este que é predominante nos turbilhóes de menor escala. Mas, 
como veremos mais adiante neste capítulo, a concorrência entre estes dois fenómenos 
(advecção e difusão) gera três escalas de movimento, a saber : energética, inercial 
e dissipativa, sendo os fenómenos mais interessantes, turbulência e intermitência, 
estão ligados a escala inercial, ou seja, a escala de transição entre a energética e a 
dissipativa. Logo, o entendimento da turbulência é a explicação da existência destas 
escalas. 

3.3.1 Escalas em Turbulência 

E uma constatação visual que em um fluxo podemos identificar estruturas de diver- 
sos tamanhos, por exemplo, nos fenómenos atmosféricos somos capazes de distinguir 
tornados, pequenos redemoinhos, deslocamentos de centros de baixa pressão. Isto 
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é possível porque porções do fluxo compartilham propriedasdes dinâmicas iguais ou 
muito parecidas e, de alguma forma, diferente das propriedades da vizinhança(fenômeno 
este ligado ao aumento abrupto do gradiente do campo de velocidade). Uma primeira 
aproximação é definir turbilhão ou redemoinho como qualquer região do fluxo, cuja 
velocidade não varia apreciavelmente. Mas agora defrontamo-nos com o problema 
de quantificar este apreciavelmente. Podemos utilizar a transformada de Fourier 
para esta análise. 

Se u{x, t) é um campo de velocidade e u{k, t) sua transformada de Fourier defin- 
imos como Turbilhão de escala menor que l = l/k a função 

u^{x,t,k)= í u{s,t)e''-''d^s. (3.30) 



|s|<fe 



Definimos, também. Turbilhão de escala l = l/k a. função 

,t,k)= J u{s,t)e''-''ds\ (3.31) 



U[X, 



\s\=k 

Definimos escala longitudinal integral ou simplesmente escala integral por 

E fácil ver, pela definição e pelo comportamento de /, que o valor da escala integral 
é determinado principalmente pela região energética do espectro. 

Definimos como região dissipativa, a região do espectro dos pequenos turbilhões, 
ou seja, dos turbilhões que são responsáveis pela dissipação da energia do fluxo 
em forma de calor. Dentro desta visão, podemos, ainda, diferenciar duas regiões: 
primeiro uma grande região, chamada de região estacionária, definida como o con- 
junto de todos os turbilhões, cuja energia E{k,t) tem dependência temporal fraca 
ou nula. Esta região, geralmente, abarca a região dissipativa e a chamada região 
inercial, que é a parte da região estacionária na qual a dissipação é pequena, ou 
seja, é dominada pelas forças inerciais (termo não-linear da NS). Esta divisão tem 
origem em algumas leis experimentais e são amplamente utilizadas em modelos 
fenomenológicos ou semi-fenomenológicos [pl|] . 
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3.3.2 Leis Experimentais 

Lei 1 (Lei dos Dois-Terços) Em um fluxo turbulento, com número de Reynolds 
muito alto, a média amostrai do quadrado do incremento da velocidade entre dois 
pontos 

E{[u,{x + r)-u,{x)f} (3.33) 

e proporcional a r^l"^ . 

Observemos que se um processo estacionário ou de incrementos estacionários é tal 
que a transformada de Fourier da média espectral do quadrado do incremento é uma 
potência 1 < n < 3, então, a média do quadrado do incremento é proporcional 
a r"~ . Ou SGJ cl, SG G cl transformada de Fourier da correlação do processo u, 
vem: 

E {ui{x ^ r)ui{x)) = J e''''Fu{k)dk. (3.34) 

Então, 

E [[ui{x + r) - Ui{x)f} = 2 J {l - é^') Fii{k)dk. (3.35) 

Fazendo Fa = C|/cp", C > O observamos que ( |3.34| ) diverge para qualquer valor 
de n, portanto processos estocásticos de segunda ordem não podem ter espectro 
em lei de potência. Não obstante, considerando somente processos de incremen- 
tos estacionários, nos quais vale ( |3.35|) , observamos que para 1 < < 3 não há 



divergência. 



E [[ui{x + r) - Ui{x)f} = Cr""-^ (3.36) 
C = 2C í (l-e''')\x\-''dx. (3.37) 



Portanto, nos casos em que vale a Lei dos Dois-terços no espaço físico, vale também 
no espaço das freqiiências a Lei dos Cinco- Terços, que vem a ser a mesma coisa. 

A região de r que respeita a Lei dos Dois-Terçosp| é denominada de região inercial. 
Encontramos esta região no centro do espectro, ocupando, geralmente, três ordens de 
grandeza, sendo seguida, à direita(número de onda maior), pela região de dissipação, 
e à esquerda(número de onda menor), pela região energética. 



tt 



Em § H há uma série de exemplos experimentais e, particularmante em uma simulação. 
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Lei 2 (Lei da Dissipação finita de Energia) 5*6, em um experimento sobre um 
fluxo turbulento, todos os parâmetros de controle são mantidos constantes, exceto a 
viscosidade, que poderá ser diminuída tanto quanto possível, a dissipação de energia 
por unidade de massa tem um comportamento consistente com a existência de um 
limite finito e positivo. 



Diferentes experimentos e simulações foram feitos para verificação desta Iei|18 



39| , sendo que a dispersão estatística varia entre 20-40%, dependendo da geometria 
e do número de Reynolds em que foram realizados. 

Portando, devemos encarar esta lei com cautela, já que não existe nenhum resul- 
tado conhecido que obrigue a sua existência. Todavia, podemos escrever a dissipação 
do fluxo como 

roo 

{e) = 2v / k'^E{k,t)dk. (3.38) 



Como a região dissipativa do espectro é a principal responsável pelo valor de 
(e), a forma detalhada de E{k,t) na região dos grandes turbilhões não influencia 
significativamente a integral acima. Isto implica que a transferência de energia 
ocorre predominantemente entre modos vizinhos, como nos mostra o modelo DIA 
de Kraichnan[^, e também que a região dissipativa é estatisticamente independente 
da região energética. Portanto, o mecanismo de transfêrencia de energia é tal que os 
pequenos turbilhões perdem toda a informação sobre as condições iniciais, gerando 
uma região cujas propriedades são independentes das condições iniciais e, além do 
mais estacionárias. Este tipo de argumentação nos leva diretamente ao modelo do 
Kolmogorov. 



3.4 As Hipóteses de Kolmogorov e suas implicações 

Na teoria de Kolmogorov consideramos um campo de velocidades incompressível, 
tal que seus incrementos sejam homogéneos, isotrópicos, e tenham segunda derivada 
contínua. Esta definição tem um aspecto generalizador, pois não exige que o campo 
de velocidades tenha sua secção espacial estacionária. O incremento do campo, 
contando também sua dependência temporal, deve ter incrementos estacionários. 
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Por outro lado, a inserção da parte temporal do campo na definição não permite 
que este varie livremente no tempo. 

3.4.1 Descrição via Incrementos 

Seja a transformação de coordenadas 

y = x' - X - u{x,t) [t' - t] (3.39) 
e o incremento do campo das velocidades 

Wi{y) ^Ui{x + y)- Ui{x) (3.40) 

pela homogeneidade, vem: 

E{wi{y)} = 0. (3.41) 

Definimos um tensor de rank 2, simétrico que nos dê a parte simétrica da correlação 
do incremento por: 

Bij(y, y') = {w,{y)w^{y') + w,{y')w^{y)} , (3.42) 
que pode ser reescrita como: 

BiAv^ y') = \ {BiÁV, y') + y') - Bij(y -y',y- y')} . (3.43) 

Pelos mesmos argumentos utilizados para a forma da covariância, no caso do campo 
de velocidades ser homogéneo e isotrópico, vem: 

E{{wiWj){y)} = yi^[Ba^{y)-B^4y)] + ô,jB^n{y) (3.44) 

y 

Bkk(y) = E {[uk(y + x) - Uk{x)f} , k ^ n, d. (3.45) 

onde o índice d denota a direção radial e o índice n a direção normal entre os pontos 
considerados. Pela homogeneidade e isotropia podemos reescrever, sem perda de 
generahdade: 

Bkk{y) = E{uk{y,Q,Of}. (3.46) 
Podemos também escrever a correlação tripla 

Bijk{a,b,c) = E {uji{a)ujj{b)uk{c)} (3.47) 
Bdddiy) = E{[ud{x + y)-Ud{x)f}. (3.48) 
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É fácil ver, pela construção dos B's, que suas relações com as funções de Kármán- 
Howarth são dadas pelas fórmulas abaixo: 

B,,{r) = 2[1 -f{r)]E{ul} (3.49) 
B^^ir) = 2[l-g{r)]E{ul} (3.50) 
Bdddir) = 6k{r)E{uiy^\ (3.51) 

Se u{x, t) descreve um fluido, então existem suas derivadas de primeira e segunda 
ordem, pois a equação de Navier-Stokes é de segunda ordem. Por conseguinte, 
existem as derivadas até terceira ordem da covariância, donde ressaltam as seguintes 
propriedades: 

Bkk{0) = O (3.52) 
drBkkiO) - O (3.53) 
d^^B^iO) = 2E{[^^uJ^{y,0,0)f}^2a^ (3.54) 
d^^BnniO) = 2£;{[9la;2(z/,0,0)f} = 26^ (3.55) 

mas, pela incompressibilidade, vem que — 2a^. 

3.4.2 As Hipóteses de Similaridade de Kolmogorov 

Hipótese 1 (Primeira Hipótese de Similciridade) Para uma turbulência homo- 
génea e isotrópica, as distribuições finitas dos incrementos da velocidade são uni- 
camente determinadas pelas quantidades: u viscosidade e por{e) dissipação média 
por unidade de massa. 

Se a família F„ depende somente de u e {e), então as correlações irão depender 
somente delas também. Todavia, as correlações dependem da distância r. Portanto, 
deve existir rj, chamado de comprimento de Kolmogorov, com dimensão de distância, 
formado por v e (s), tal que: Bkkir) = Pkk{r / rj)!!^ , k = d,n, onde Pkk é uma função 
universal c u c uma constante com dimensão de velocidade que só pode depender 
de v e (e). Fazendo uma análise dimensioal simples, chegamos a: 



iyy{e)\^^^ (3.56) 
u = (3.57) 
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A turbulência é isotrópica, logo B^kij) = Bkk{—r), e por (|3.4.1| ) a função é par. 
Retendo o primeiro termo da expansão em série de potência e utilizando novamente 



Bddir) 
Bddir) 



r 



,k = d,n 



V 



(3.58) 

(3.59) 
(3.60) 



Hipótese 2 (Segunda Hipótese de Similaridade) Para uma turbulência homo- 
génea e isotrópica, quando \y \ ,\y' \ , \y' — y\ rj, as ditribuições finitas do incre- 
mento da velocidade são unicamente determinadas por (e). 

Tomando r suficientemente grande, as funções universais Pkk e Pddd devem ter 
uma forma tal, que as funções de correlação Bkk, Bddd não apresentem dependência 
em u. Portanto, pela Primeira Hipótese de Similaridade e pela ( |3.4.2| ), teremos: 



2/3 



Pkk -'^[-j ■ (3-61) 
Utilizando ( p.4.1|) e fazendo um procedimento análogo para a correlação tripla: 

B,, ^ C((e)r?/' (3.62) 
''-C{{e)rf^ (3.63) 



B, 



M = C{{e)r) 

4 
-( 

3 

C (e) r, 



R — 



ddd 



(3.64) 



onde C, C são constantes adimensionais a serem determinadas. 



3.4.3 As Hipóteses de Kolmogorov e a Relação de Kármán-Howarth 



Substituindo ( |3.4.1| ) em ( p.22| ), obtemos a RKH para as correlações escalares de 
Kolmogorov. 



3^^ 
-^dt 



2\\u\\-B. 



dd 



^^Bddd + ^^ 



( d"^ A d \ 

Qu —17 Bdd + -TrBdd , 
V ar"^ r or I 



lembrando que a turbulência é estacionária e (e) — ^ ^ " 



(3.65) 

1^ llMii , reescrevemos a equação 



acima como: 



4{e) 



d 4 
dr r 



d 



Qi^T^Bdd — Bddd , 
ar / 



(3.66) 
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integrando e usando as condições iniciais vem: 

4 d 

-{e)r = 6u—Bdd - B^d- (3.67) 

5 ar 

Observemos que, se tomarmos r da ordem do comprimento de Kolmogorov, 



esperamos que Bddd ^ e por ( |3.67| ) temos que §^Bdd\r=o ~ ji^r, ou seja B, 



dd 



jl^r^. Tomando r muito maior que o comprimento de Kolmogorov, esperamos pelas 
Hipóteses de Similaridade, i.e. ( p.4.2|) , que ^Bdd ^ Bddd, portanto Bddd = | (^)^- 



3.5 Objeções à Universalidade de Kolmogorov 

Landau foi o primeiro que suspeitou da universalidade das constantes, baseado no 
seguinte argumento: 

Suponhamos, a princípio, que o valor instantâneo de seja expresso como 

uma função universal da dissipação e(t). Quando tomarmos a média desta ex- 
pressão, a variação de e sobre tempos da ordem do período dos turbilhões grandes(tempo 
carcterístico de advecção) influencia o valor da média. Todavia, esta variação é di- 
ferente para cada fluxo, pois a região energética é permanente e caracterizada pelas 
condições de contorno do problema. Logo, o resultado não pode ser universal. 

Este argumento também é utilizado por Kraichnan||36||, só que de uma maneira 
mais detalhada e quantitativa. Cuja linha mestra é: 

A variação temporal de e (t) é determinada pela convecção da escala dissipativa, 
pelas escalas maiores (pertencentes à região inercial) e por suas distorções internas, 
sendo que a própria convecção pode causar distorções. E por seu turno, as dis- 
torções internas e o período de convecção são características de cada fluxo, assim as 
constantes para as correlações não têm um caráter universal. Então, se e independe 
dos detalhes do fluxo, a região inercial deve ser estatisticamente independente da es- 
cala de dissipação. Ou seja, o mecanismo de transferência de energia dos turbilhões 
deve apagar detalhes das informações contidas nos turbilhões maiores e ter um 



caráter local. 

Todavia, Kraichnan mostrou que, se considerarmos inicialmente a escala dissipa- 
tiva estatisticamente independente da região inercial e deixarmos que a única forma 
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de interação entre elas seja o termo não-linear estas escalas passarão a ser esta- 
tisticamente dependentes e a ter momentos depententes do tempo para instantes 
diferentes. 

Para ilustrarmos os parágrafos anteriores, calculemos o exemplo de Kraichnan : 

Exemplo 3.1 (Kraichnan) Tomemos um campo de velocidade u + v tal que: v 
seja constante no espaço e no tempo; tenha uma distribuição gaussiana e isotrópica; 
u seja variável no espaço, muito menor que v, e em í = O tenha ditribuição gaus- 
siana, homogénea, isotrópica e independente de v. Escrevendo a equação de Navier- 
Stokes sem o termo viscoso, oomo o termo viscoso é responsável pela dissipação 
e o fenómeno em questão é a interação, este somente atrapalharia. Assim, cada 
componente de Fourier da velocidade obedecerá à seguinte equação: 

d 

—u{k,t) = i{k ■ v)u{k,t) 

u{k,t) = uik,0)e-'^''''^\ (3.68) 
Calculando a correlação temporal 

E{u{k,t)u*{k,s)} 



T{k;t,s) 



[E {u{k, s)u*{k, s)} E {u{k, t)u*{k, t)}f^^ 



onde Vo é a variância de v. Mas, por hipótese, u e v são independentes no instante 
inicial, logo os valores simultâneos dos dois campos são estatisticamente indepen- 
dentes em qualquer instante posterior. Decorre de ^3. 6dJ que 



E {v,uj{k, t)u*Jk, s)} = -ivlh {t - s) e-H>^'^'-^?E {u,{k, OX(A;, 0)} (3.69) 

Poderíamos generalizar o resultado acima e mostrar que momentos mais gerais so- 
mente se anulam para tempos iguais, tal qual o exemplo acima. 

3.6 Observações 

Muito embora a Lei dos Dois- Terços se aplique a uma grande variedade de fluxos, 
como nos mostrou Kolmogorov, esta não parece estar completamente de acordo 



com os resultados experimentais ||18[| , já que estes parecem indicar a existência de 
termos de correção, ou seja, i?„„ ~ ^2/3+ç_ yig^^ ainda não temos dados, com pre- 
cisão suficiente, para saber se C é causado pela inaplicabilidade das Hipótese de Kol- 



mogorov ou pelas técnicas intrusivas utilizadas na maior parte dos experimentos p9[. 
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Poderíamos pensar em realizar simulações para termos comprovações, mas para este 
fim seria necessário um tempo de máquina impensável. 



Existem, hoje, evidências da existência de estruturas coerentes |TÕ|, ^ nas soluções 
da equação de Navier-Stokes, como veremos mais adiante para o caso da equação de 
Burgers, estas estruturas são responsáveis pelas deformações da distribuição de pro- 
babilidade. Assim, ( pode ser o resultado da existência destas estruturas na região 
inercial. Contudo, a quantificação desta influência não é trivial devido a interação 
entre turbilhões de diversos tamanhos. Mas as simulações vêm avançando muito nos 
últimos anos, e no futuro próximo teremos uma resposta finalQ para esta questão 
secular. 



^Começará, no próximo ano, um esforço conjunto entre o Observatoire de la Côte D'Azur(Prof. 
Frisch) e o Laboratório de Los Alamos(Prof. Kraichnan) para executar uma simulação cuja malha 
será da ordem de 2048"^ ^ 10^^ pontos. 
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Formulação Funcional 



A análise da Turbulência via funcionais, mais precisamente, por meio do funcional 
característico da distribuição de probabilidade do campo de velocidades, começou a 
ser aventado por Kolmogorov. Todavia, sua primeira formulação sistemática deve- 
se a Hopf|2^, que introduziu uma equação integro-diferencial para o funcional 
característico. Mas esta abordagem é incompleta, pois somente os momentos espa- 
ciais, ou seja, momentos calculados no mesmo instante de tempo, são passíveis de 
serem calculados. Lewis e Kraiclinan[[41| generalizaram a abordagem de Hopf de 
tal forma que os momentos espaço-temporais podem ser calculados. Como o pro- 
cesso dedutivo das duas formulações é análogo, utilizaremos como base o artigo de 
Lewis e Kraiclinan. Cabe ressaltar aqui que, embora a equação de Hopf não nos 
forneça uma descrição completa do fenómeno estatístico, ela está em melhor posição 
para ser comparada com os resultados canónicos do modelo de Kolmogorov já que, 
quando tomamos o limite de tempos iguais para os momentos espaço-temporais, 
estes tendem, forçosamente, para zero se a separação espacial permanece finita. 



4.1 Funcional Característico 

Definimos o valor médio de um funcional F sobre o espaço de soluções da equação 
de Navier-Stokes para uma determinada condição de contorno, como: 

(F) = j dP{u)F{u), (4.1) 
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onde P (ti) é a medida de probabilidade neste espaçoQ. 

Agora seja y = {yi, 112,113) = y{x,t) um campo real arbitrário, que se anula no 
infinito espacial, em particular, podemos tomar y como tendo suporte compacto. 
Definimos um funcional de y, que é Funcional Característico de P, como sendo: 

F[y] = (exp2 {/^^ dt J^^ d^xyi{x,t)ui{x,t)^'j = (e^^^'^)) . (4.2) 

Calculando a n-ésima derivada funcional parcial de F para y = O teremos a seguinte 
igualdade: 

^ 7 7 \y=0= {UaÁXutl)---Ua„{Xn,tn)) . (4.3) 

Por outro lado, caso desejemos obter a função característica da densidade de proba- 
bilidade a n-pontos basta tomarmos 

n 

yi{x,t) = hiQ)ôia^5{t - tj)5^{x - xj), (4.4) 
i=i 

onde em ^i(j), o primeiro índice indica a componente do campo e o segundo o ponto 
referido. Substituindo na definição do Funcional Característico temos 



n 



F[y] = ( expi \ K,(j)Ua,{xj,tj) ) ) = p(/i(i), /i(„)), (4.5) 



logo, a densidade de probabilidade é 



p{Uai{Xi,ti), ...,Ua„{Xn,tn)} = (4.6) 

» n 

(27r)""y d/i(i)...c//i(„)p(/i(i),..., )exp{-i^M„^./iQ)}. (4.7) 



4.2 Equação Para o Funcional Característico 

Para obtermos uma equação diferencial para F, basta observarmos que 

5F[y] 



Sya{xi,ti] 



í(M«(xi,ti)e*(^'")) (4.^ 



* Observe que este funcional é essencialmente diferente do encontrado em Mecânica Quântica, 
onde as integrais para as probabilidades de transição podem ser em alguns casos discretizadas e 
depois tomado o limite do contínuo e assim calculadas diretamente. 
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Sya{xi,ti)ôyi3{x2,t2) 

_d d SF[y] 

dxp dxp Sya{x,t) 

d S^^Flv] 



d d 



dxp dxp 



d 



-{e''^y'^'>up{x,t)-^^ua{x,t) 



(4.9) 
(4.10) 
(4.11) 



dxp Sya{x,t)5yp{x,t) 

Para obtermos a última equação, utilizamos a condição de incompressibilidade. As- 
sim, também teremos 



d ÔF[y] 



i ( e 



i{y,u) 



^Ua(x,t)). 



(4.12) 



dtôya{x,t) " y dt' 

Substituindo, no lado direito da equação (4.12), a equação de Navier-Stokes com 
uma força determinística e utilizando as relações acima vem: 



d 5F[y] _ . d S^F[y] 



dt5ya{x,t) dxpSyc,{x,t)õyi3{x,t) 
d d SF[y] 1 9 ^ 



(4.13) 



dxf3 dxfs 5ya(x,t) ' '''' dxa 

onde n = i (p{x,t)e^^y'^'''^. Para eliminarmos o campo dc pressão, ao invés de 
recorrermos à utilização de operadores não locais, preferimos introduzir um campo 
de teste l{x,t) — {li,l2,h), tal que este se anule suficientemente rápido no infinito 
espacial e seja solenoidal. Fazendo uma integração por partes vem : 



J \âxi j 



- I d^xU í^k 

\ OXi , 



0. 



(4.14) 



Assim, integrando o campo l conjuntamente com a equação para F, obtemos 

5^F[y\ 



j dtd^xla 



d 5F[y] _ . d 



dt5ya(x,t) dxp5ya(x,t)5yf3{x,t) 



(4.15) 



j dtd^xla 



dxpdxp Sya{x,t) 

Lembrando que esta equação precisa ser satisfeita para qualquer campo de teste 
solenoidal e pode ser tomada como substituta da equação de Navier-Stokes, com a 
condição inicial reescrita como: 



Po[y] = (exp<^i / d^xyi{x,0)ui{x,0) 



(4.16) 
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A condição de incompressibilidade pode ser escrita como: 

âXa Òya{x,t) 

que é uma conseqiiência direta da definição de F. Cabe lembrar também que F tem 
que ser normalizado e a probabilidade deve ser não negativa, implicando as seguintes 
relações 

F[0]-1 F*[y]^F[-y] F[y] < 1 (4.18) 



E 



i{yj,u) 



J2J2^>jF[yj-yi\>Q- 

1=1 j=i 



(4.19) 



Esta última propriedade é chamada de definição positiva do funcional. Podemos es- 
tabelecer, também, uma equação para o funcional característico utilizando a trans- 
formada de Fourier dos campos determinísticos, ou seja, efetuando as seguintes 
transformações : 



Zi {k, t)^ j d^xyi {x, t) e 
yi{x,t) = {"^T^y^ j d^kzi{k,t) e 



(4.20) 
(4.21) 



rrii (/c, t) = I dPxli (x, t) e '''^^^ 



k (x, t) = (27r) ^ j (fkrrii {k, t) 



Qi {k,t) = (27r)-'y d^xfi{x,t) e 
fi (x, t)^ j d^kgi {k, t) 



Como y e f são reais, segue da transformação de Fourier que 



(4.22) 
(4.23) 

(4.24) 
(4.25) 



z* {k,t) ^ Zi{-k,t) e g* {k,t) ^ gi{-k,t) , 



(4.26) 



donde vem, também, que a incompressibilidade de li fica representada por 



kiruj — 0. 



(4.27) 
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Esta transformação no espaço das funções implica que 



SF[y] f^3^^ ÔF[z] 5zpik,t) 
6ya{x,t) J Szp{k,t) õya{x,t) 

Portanto, para a n-ésima derivada teremos 



d^k 



5F[z] 
Szp{k,t) 



(4.28) 



Ô^F[y] 



(4.29) 



d^ki...£kr, 



Ô^^Flz] 



Zi, ti) . ..ÔZct^ (^nj tn) 

Substituindo na ( [4.15|) veni|: 



SyQ-,{Xi,ti)...Sya„{Xn,tn) 

exp |-i xi^ki^ + ... + Xn^kn^ I . (4.30) 



y dtd^krrii (/c, t) 



d_ 
di 



+ vk^ 



ÔF[z] 
ózi{k,t) 



- íF [z] 



kikj 
5 a r^e- 



dtd kid k2mi (fci + k2,t) ^1+^2 



9j {k,t) 
Õ^F[z] 



(4.31) 



ôzi{ki,t)ôzj{k2,t)' 

A equação acima nos sugere que podemos generalizar nossa abordagem um pouco 
mais tratando a força como sendo um processo estocástico. 

Suponhamos que no instante inicial a velocidade e a força, sejam estatisticamente 
independentes e dados. Logo, ao invés do funcional F[y], consideraremos um novo 
funcional que descreva a probabilidade conjunta da força e da velocidade. Assim, 
definimos 

F [y, q]=JdP {u, f) exp z 11 dtd^x {y^m + g./^) | . (4.32) 
Agora, se tomamos 



yi (x, t) = 6(t-0) Ti (x) 
qt (x, t) = ô(t -0) Si (x) . 



(4.33) 
(4.34) 



substituindo na definição de F e utilizando a condição de independência, vem: 



F [y, ç] = ( exp i / d^x [r^ (x) Ui (x, 0) + Si (x) fi (x, 0)] 



dP (u) exp i / d^xri (x) u (x, 0) 



X 



dP (/) exp i / d^xsi (x) fi (x, 0) 



TPara obtermos esta equação, utilizamos o operador de projeção 
que o campo de velocidades responde somente à parte solenoidal da força 



ç kikj 



X [r] ^ [s] , 

para ressaltarmos 
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onde X M [s] são respectivamente os funcionais característicos para a velocidade 
e para a força no instante inicial. 

Fica fácil ver que a modificação que incide sobre a ( |3.13| ) é simplesmente 

5F[y,q] 



ifaF[y,q] 



levando-nos a reescrevê-la como: 

d 6F[y,q] 



ôqa {x,ty 



d Ô^F[y,q] 



dtõya{x,t) dxi3Syc,{x,t)õy/3{x,t)_ 
dtd^xl 



^, d d ÕF[y,q] ^ ÕF[y,q] 



dxfs dxfs õya{x, t) ôqa {x, t) 
sendo que a condição de incompressibilidade não muda formalmente. 



(4.35) 



(4.36) 



4.3 Solução para Número de Reynolds Nulo 

Observar o comportamento da equação para o número de Reynolds nulo significa que 
estamos olhando para uma equação onde os efeitos de transporte de energia entre 
os turbilhões são desprezíveis, ou seja, toda a energia recebida da força externa é 
dissipada através do termo viscoso. Em outras palavras, esperamos que esta solução 
descreva o comportamento da região dissipativa do espectro. 

Comecemos pelo caso mais simples em que a força externa é determinística. 
Obviamente o termo que deve ser eliminado é o termo na segunda derivada funcional 
de F, que na representação espectral é o termo essencialmente não local. Portanto, 
a equação toma a seguinte forma: 

5F[z\ 



dtd krrii (k, t) 



9 ;2 



- íF [z] 



kikj 



gj{k,t)\=0. (4.37) 



õzi{k,t) 

Como a integral acima deve ser nula para todo campo m, necessariamente a 
expressão entre parênteses deve ser nula, ou seja, 

ÔF[z] 



9 ;2 



- íF [z] 



kikj 



9jik,t) = 0. 



ôzi{k,t) 

Suponhamos que o termo de força seja nulo. Logo, a equação fica reduzida a 

ÕF[z[ 



(4.38) 



d ,2 



õzi{k,t) 



0. 



(4.39) 
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Como a equação acima nao tem derivadas funcionais em z, devemos ter: 



ÔF[z] 



Ri [z] e 



-vkH 



(4.40) 



5zi{k,t) 

Como neste ponto o funcional Ri [z\ é arbitrário, tomemos Ri [z\ = Zi{k,t). Mas, 
por ser F um funcional característico de fluxo incompressível, devemos ter: 



kiZi — 0. 



(4.41) 



Agora, como queremos tomar z como um campo arbitrário, basta inserirmos o ope- 
rador de projeção solenoidal diante do campo z na nossa solução. E para descrever- 
mos uma condição inicial arbitrária, já que Ri [z] deve ser essencialmente a derivada 
funcional do funcional que é solução, escrevemos a solução como: 



F[z] = A 



dte 



Zj{k,t) 



(4.42) 



Para resolver ( |4.38 ), podemos tentar multiplicar o funcional, que já sabemos ser 
solução da equação não forçada, por um funcional que ao ser derivado cancele o 
termo forçado. Ou seja, tentemos substituir 



F = A- 



dte 



5.. _ -^im 
k^ 



Zj{k, t) > exp 



i I dtd^khjZj 



(4.43) 



na equação forçada, levando a seguinte equação: 

d 



dt 



hi {k,t) 



kikj 



9j {k,t) 



(4.44) 



equação que tem como solução 

hi {k, t) = e 



'^^'^ ' dre 



„ kj k A 

A. . LJ. 



9j {k, t) 



(4.45) 



Consideremos, neste momento, o caso em que a força é aleatória e o número de 
Reynolds nulo. A ([4.36| ) fica 



d 5F[y,q] d ÔF[y,q] _ 5F[y,q] 



dt6ya{x,t) dx(3dxp6ya{x,t) 6qa{x,t)' 
Transformando para a representação espectral, onde 



(4.46) 



g,ik,t) = {2n)-' I d:'xq, {x,t)e''^^^ 



(4.47) 
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teremos 



9 ,2 



SF[y,q] 

ÔZa{k, t) 



a(3 



(4.48) 



Por analogia com ( 4.43| ) e lembrando que a condição inicial é que o campo de força 
seja estatisticamente independente do campo de velocidade, teremos como solução 



X y dte 
X \l/ < I dte 



vk^t 



kikj 



Zj{k,t)^ X 
Zj{k,t) +g{k,t) 



(4.49) 
(4.50) 



4.4 Abordagens Ingénuas de Solução 



Possuímos uma equação diferencial para o funcional característico que, à primeira 
vista, nos parece muito aprazível, pois partimos de uma equação não-linear e in- 
omogênea e chegamos a uma equação linear e homogênea( |4.36| ) . 

Todavia, não existe um método geral para resolver equações funcionais integro- 
diferenciais, nem tampouco teoremas de existência e unicidade das soluções para 
o caso geral]]. Outra dificuldade é que equações deste tipo não são de fácil imple- 
mentação computacional. 

Portanto, nesta secção, tentaremos mostrar algumas abordagens já realizadas e 
onde aparecem as suas dificuldades. 

Reescrevendo a equação para o funcional característico como 



^^DyF[y, q] - iVDyDyF[y, q]^ = 
il,iyADyF[y,q] + D,F [y,q]) 



(4.51) 



onde (a, b) = J dtd?xaibi é o produto interno do espaço das soluções da equação de 



Navier-Stokes, D 



^ - ^ e D = ^ 

' Sy(x,t) ^ ^1 Sq{x,t) 



sao "gradientes funcionais" em relação 



aos campos y {x, t) e q (x, t), respectivamente, sendo que os outros símbolos têm seu 
significado usual. 



■'■Recentemente tem havido um grande esforço no estudo de equações diferenciais funcionais em 
espaços de Banach. Mas nenhum resultado significativo é de nosso conhecimento (wde; Approx- 
imation Methods in Banach Spaces, Mesure Theory in Infinite Dimension Spaces MIA, Kluwer 

)■ 
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Em uma primeira vista, poderíamos pensar em expandir F como uma série de 
funcionais à la Taylor, ou seja, 



oo oo 



F[y,g]=F(0'°) + 5:5:F("'-), (4.52) 



n=l m=l 

onde 



p{n,m) ^ j dt^d^Xi... J dtnd^Xn j dtn+ld^Xn+l-. j dt^+md^Xn+m X (4.53) 



àya^ (Xi, ti) ...Ôya„ {Xn, tn) (x^+i, í„+i) ...í 



(4.54) 



X yai {Xi,ti) ■■■yan {Xn,tn) X Ça^+i (^^n+l) ^n+l) ■■■Ça„+rn {Xn+mjtn+m) ■ (4.55) 

Assim, é evidente que 

= F [O, 0] = 1. (4.56) 

Mas desta forma, ganhamos o problema de convergência da série, já que os 
momentos não necessariamente se anulam quando, pelo menos dois pontos estão 
separados por uma distância infinita. 

Se, por exemplo, ignoramos o termo de força de ( [4.51| ), fazemos l = y, utilizamos 
a expansão proposta em (|4.52| ) e agrupamos os polinómios de mesma ordem em y, 



obtemos 

'd 



D,F^"'] = {v, !VB,D,Fl"+'l) , (4.57) 



que é exatamente a equação para o momento da n-ésima ordem, sendo que este, 
como já havíamos visto, envolve o cálculo do momento de ordem subseqiiente. 

Portanto, a única vantagem, neste caso, é estética, já que deixa a equação de uma 
forma mais compacta. O problema de convergência da série pode ser contornado, 
utilizando uma expansão no logaritmo de F, ou seja, uma expansão nos cumulantes, 
garantindo assim, a convergência das integrais, mas obviamente não resolvendo o 
problema de fechamento. 
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Poderíamos tentar utilizar uma expansão do tipo Gram-Charlier|, ou seja, 



F = e 



i?(l)+i?(2) 



1 + g 



n=3 



(4.58) 



Mas esta também nos leva a equações que precisam de fechamento, muito embora 
a determinação destes F*^") seja equivalente à determinação de um número infinito 
de momentos. 

Neste ponto, poderíamos aventar a possibilidade de truncar a expansão em uma 
ordem específica, para conseguirmos uma primeira aproximação do problema. Mas, 

k oo 



n=k+l 



n=l 



< 1 









oo 


\F\ < 


1 + ^ 


+ 






n=l 




n=k+l 



(4.59) 
(4.60) 



Se supusermos que Y^'^=k+i F^^^ = O, então, existirá pelo menos um y, a saber 
y = O, tal que 1 + Y.n=i -^''"^ > ^- Logo, a série truncada não pode ser Funcional 
Característico de nenhum processo estocástico. Ou seja, só é válida a expansão 
de Taylor com resto, não nos fornecendo, portanto, um método de aproximações 
sucessivas. 

Em outras palavras, mesmo que não saibamos exatamente todos os momentos, 
cumulantes ou desvios gaussianos(dependendo da expansão utilizada), temos que ter 
uma estimativa para eles, se quisermos utilizar este método. Contudo, estimativas 
normalmente virão de uma hipótese ou modelo, e não diretamente da equação para 
o funcional. 

Poderíamos, ainda, tentar fazer uma expansão da solução com base no número 
de Reynolds mas este procedimento é desaconselhado pelos seguintes aspectos: (1) 
estamos interessados em Turbulência Desenvolvida, logo devemos tomar o limite 
quando o número de Reynolds vai para o infinito. Este fato obriga-nos a fazer 
ressomas de séries perturbativas, cuja convergência não é bem estabelecida, podendo 



§A conveniência desta expansão advém do seguinte fato: se ignoramos os termos da segunda 
derivada funcional, ou seja, Número de Reynolds nulo, e de força, a solução será estacionária, o 
que significa que, se tivermos um campo gaussiano, este permanecerá como tal para todo tempo 
posterior [El]. 
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levar-nos a resultados errados em caso de convergência ||1 7|| ; (2) a invariância por 
transformação aleatória de Galileu é quebrada para qualquer aproximação, pois 
a expansão é feita sobre o termo não linear da equação, sendo que a invariância 
é conseguida através da anulação do mesmo com o termo de derivada temporal; 
(3) todos os resultados conseguidos por métodos diagramáticos foram conseguidos 
também por métodos mais simples. 

Mesmo com estas dificuldades, alguns físicos vêm utilizando métodos diagramáticos 



com um relativo sucesso( ver[^] e suas referências), muito embora nos pareça im- 
provável que este método consiga obter resultados que não foram obtidos por outras 
técnicas menos abstratas. 



Capítulo 5 



A Equação de Burgers 



Neste capítulo trataremos da equação de Burgers unidimensional, que é uma redução 
para uma dimensão da equação de Navier-Stokes, sem o termo de pressão. Esta 
equação, embora introduzida por J.M. Burgers[^ como um modelo para o estudo 
da turbulência, não apresenta caos, mesmo na versão forçada. Por outro lado, ela 
surge naturalmente em estado sólido no estudo de polímeros diretos e crescimento 
de superfície, neste caso sendo mais conhecida como equação KPZ [p^] . 

O interesse neste modelo simples advém do fato de ser integrável (via trans- 
formação de Hopf-Cole[^) e termos muitas informações sobre as propriedades das 
soluções, no caso determinístico [0, |T3|, 0, ^ e uma formulação matemática pre- 



cisa para o caso estocástico ||T9l. Isso faz com que este modelo seja um laboratório 
ideal para o teste de técnicas que posteriormente devem ser aplicadas à equação de 
Navier-Stokes. 

Este capítulo tentará dar uma visão panorâmica sobre a equação determinística, 
e mostrar os principais resultados conhecidos para o caso aleatório. 

5.1 A Equação Determinística 

A equação de Burgers é 

d d 

= -M— M + Z/— M, (5.1) 

também podendo ser escrita na forma dita conservativa 



d_ _d_ 

dt dx 



(5.2) 



54 



Capítulo 5. A Equação de Burgers 



55 



quando adquire uma forma mais palatável, se introduzimos a transformação de 
Hopf-Cole, 



d 

u = ——ip, ip = 2vhiw, 
ox 

levando-nos, a menos de uma fase arbitrária no tempo, à equação de calor 



d 



(5.3) 



(5.4) 



Se tomarmos como domínio toda a reta real e tempos estritamente positivos, a 
solução é 

(5.5) 



Escrevendo a condição inicial como 



w (x, t) = / e wq [a) da. 



Wçj (x) = lirn w (x, t) , = e-^^^'' 



vem, da transformação de Hopf-Cole, que 

Uq (x) = - 



ax 



u (x, t) 



2t 



Jdie 



(5.6) 



(5.7) 



(5.8) 



Importante observar que, para esta solução existir, devemos necessariamente ter 

/o^áC«o(C) 



lim ^„ 
5.2 Aproximação Heurística 



0. 



(5.9) 



Para viscosidade nula, a equação de Burgers, dada a existência do Teorema de 



Hopf|^, pode ser escrita como: 

d 



dt 



, , ^ d d d 



(5.10) 



ou seja, as partículas lagrangeanas conservam suas velocidades. Assim, a solução 
desta equação pode ser escrita como: 



u (X (a, í) , t) = Uq (a) 

X(a, t) = a + íuo (a) 



(5.11) 
(5.12) 
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onde utilizamos letras maiúsculas para coordenadas eulerianas,i.e. u {X (a, t) ,t) e 
letras minúsculas latinas para coordenadas lagrangeanas, i.e. a. Todavia, esta é uma 
solução implícita da equação de Burgers, chamada de lagrangcana, que permanece 
única, inversível e regular, até o instante que o jacobiano da transformação euler- 
lagrangeana permanece inversível, ou seja, até que 

d 

J{a,t)^—X{a,t)^l + túo{a) (5.13) 

seja positivo. O instante em que isto acontece é denotado por í*, e o chamaremos 
de tempo de pré-choque , sendo definido por: 

t* = ■ V . .y (5.14) 

-mm^ [uo {a)\ 

Notemos que no caso em que Uç, (a) — ca c > O, o tempo de pré-choque é negativo, 
ou seja, se no instante inicial não existe nenhuma partícula no fluido capaz de 
ultrapassar a sua vizinha, o fluxo permanecerá sempre numa fase regular (onde todo 
o fluxo forma um único choque). Isto pode ser melhor visto sc substituímos esta 
solução na equação para a coordenada euleriana, i.e. X (a, t) — {1 + te) a. Como 
soluções não decrescentes, no conjunto das soluções da equação de Burgers, são 
raras, não precisamos concentrar nossa atenção sobre este caso. Concentraremos 
nossa atenção no conjunto de soluções cujos gradientes sejam não constantemente 
positivos. 

Seja a* a posição lagrangeana tal que í* = l/[— Wo(a*)], portanto a derivada 
espacial da velocidade euleriana em um instante t, e no ponto X^ (a*, t) será 



d du da . . 1 /„ ^ 

que é inflnito no tempo de pré-choque. Por isso, chamamos esta singularidade em a* 
de pré-choque. Fazendo uma expansão de Taylor ao redor do pré-choque, tomando 
«o (o*) = 0,a^ = O, condição esta que pode ser conseguida por transformação de 
coordenadas e transformação de Galileu, vem: 

1 b 

Uo{a) = a + -a^ + 0(a^). (5.16) 
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Como o gradiente da velocidade é mínimo global em a*, devemos ter a segunda 
derivada nula e a terceira positiva, ou seja 6 > 0. A expansão para a posição 
euleriana fica 

X(a,í) = (l--)a+^a3 + 0(a^). (5.17) 

Se desejamos inverter esta transformação, teremos dois comportamentos diferentes 
na vizinhança do pré-choque, dependento do tempo, ou seja, 

1 



a (X, t) 



a {X, í*) 



X + (x^) , t<u 



(5.18) 
(5.19) 



implicando, para a velocidade, na seguinte expansão: 

[t* t) 

1 /3 

u{X,t) = -^(^^) t^U 



(5.20) 
(5.21) 



Assim, podemos visualizar o pré-choque como uma mudança local de um compor- 
tamento linear decrescente, í < í*, para um comportamento tipo raiz cúbica, t — t*, 
com tangente vertical. 

Podemos ainda analisar o comportamento da enstropia que é definida como: 



fl{t) 



1 



d_ 
dx 



u {x, t) 



dx 



I 



Uo 



1 + túo 



[1 + túo) da, 



(5.22) 



obtendo 



Q (í) = oc (í, -t) í t U 



n(t) = oc I \x 

leio 



-4/3 



(5.23) 
(5.24) 



Estando o pré-choque caracterizado, vejamos como ele se forma. Para tal, con- 
sideremos que o campo inicial é analítico, tal que possamos fazer a extensão da 
solução implícita ao domínio complexo, denotando U {Z, t) a velocidade euleriana 
complexa, A a coordenada lagrangeana. Assim, se a aplicação complexa 



Z^A + tUo {A) 



(5.25) 
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deixa de ser inversível em um ponto A* (t), cuja imagem pela aplicação é (t), 
então a expansão em torno deste ponto será 



Z = oc Z^t) + {A - A4t)y 



dZ 



dA 



oc \A-A,{t)\ocJ\Z-Z,{t)\ 



e, portanto, o comportamento local do gradiente da velocidade será 



dU 



dZ 



dU 



dA 



dZ 



dA 



(xl/J\Z - Z4t)\ 



U-U,= oc J\Z-Z4t)\ 



(5.26) 
(5.27) 

(5.28) 
(5.29) 



ou seja, a existência de singularidades complexas no campo inicial que se aproximam 
da reta real, conforme t se aproxima de t^,, dá origem aos pré-choques. 

Para entendermos globalmente a solução, comecemos por reescrever ( |5.8|) como: 

^ t 



u (x, t) 



onde 



(x, t, a) 



h / adãuQ (ã) 

2t Jo 



(5.30) 



(5.31) 



Em uma solução para u tendendo a zero, os únicos pontos a que contribuem para 
( ^.30|) , para x,t fixos, são os pontos nos quais (p{x,t,a) atinge um mínimo global. 
Ou, de uma maneira mais formal, definimos 



ip (x, t) = min (x, t, a) 



que pode ser reescrita como 



X 



onde 



'^{x,t) = ^ + Lt (a(x,t),) 



r xd '] 

Lt{a{x,t)) = min|u;(x, a) ^| 

uj{a,t) = 



da 



uo [a) + 



(5.32) 

(5.33) 

(5.34) 
(5.35) 



Definimos uma função (a), chamada de envelope convexo, tal que goze da pro- 
priedade de ser a maior função convexa, tal que < u, podendo ser descrita 



Capítulo 5. A Equação de Burgers 



59 



geometricamente, por intermédio da seguinte construção: fixando x et, construímos 
uma reta no plano {uj, a), que seja descrita pela equação u = + c. Para cada x 
podemos encontrar cq {x, t) = Cq tal que, para todo c < Cq a reta uj = ja + c não 
intercepte o gráfico de cu, enquanto que para todo c > Cq a interceptação acontece. 
Mas para c — cç, a, reta u — |a + Cq é tangente ao gráfico de cu, em pelo menos 
um ponto. Portanto, ip {x, t) é, em geral, um conjunto dos pontos a, onde a reta 
cu — jã + Cq é tangente ao gráfico de cu. Definimos, então, 

a-{x,t) = min {a; a E {x,t)} (5.36) 
a+{x,t) = max{a;aEip{x,t)}, (5.37) 

ou seja, estes são o menor e o maior valor da aplicação plurívuca ip {x, t). 
Finalmente, definimos o envelope convexo como: 

\ uj{a,t), a.{x,t) = a+{x,t) 
Coj [a) = < (5.38) 
I + Co, a_ {x, t) < a+ {x, t) . 

Consideremos, agora, sua derivada Ft (a) = (a), que em geral será uma função 
não-dccrescente, e sua inversa Ff-^ (x), que geralmente não é uma aplicação bem 
definida em razão dos pontos x, para os quais existe um conjunto {a; Ft (a) — x}. 
Nos pontos em que Ff^ (x) for descontínua, vale a pena considerar Ff^ (x) como uma 
curva contínua sobre o plano , que tenha segmentos verticais para estes valores de x, 
já que posteriormente estaremos interessados em uma interpretação mais pictórica. 
Agora, podemos formular o seguinte teorema: 

Teorema 5.1 (Hopf) Sejam x,t tal que %l){x,t) consista de um ponto a{x,t) — 
Ff^ (x) . Então lim^_>o u{x,t) — Uq (x) existe e 

X — Fr^ ix) . ^, 

u {x, t) = (5.39) 

t 

se F~^ {x) é um intervalo de comprimento positivo, então existem os limites laterais 

ií_ {x, t) = lim tio (x + x) ^ - — "~ (5.40) 
í— >o~ t 

u+{x,t) = fim uo{x + x) ~ ""^ . (5.41) 

í— +0+ t 
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Ao invés de fazermos a demonstração matemática deste teorema, vamos dar uma 
visão mais "física" dele. 

Como a (x, t) é a coordenada da partícula, onde (x, t, a) atinge um mínimo 
global para um dado é fácil ver que a{x,t) é a coordenada lagrangeana da 

qual emana a partícula do fluido, que estará na coordenada (x,t) euleriana. 

Assim, podemos interpretar o resultado de Hopf e os resultados heurísticos da 
seguinte forma: o choque entre partículas lagrangeanas do fluxo, originadas da 
variação espacial da velocidade, gera a coalescência entre estas partículas, refletindo- 
se na perda de inversibilidade da transformação de coordenadas eulerianas para la- 
grangeanas, formação dos pré-choques(flg: ^.21) , implicando, dessa forma, na perda 
de continuidade do campo nestes pontos, sendo que estes pontos de descontinuidade 
são chamados de choques(flg: f)-l\)f\ Assim, o perfll de velocidades do campo eule- 
riano é tipo serrote, sendo esta fase chamada de fase dissipativa. Como ilustração, 
consideremos o caso abaixo. 

Suponha agora, que a solução u {x, t) apresente somente um choque, localizado 
em X (t), e que 




Figura 5.1: Formação do choque para condição inicial periódica. 



* estas figuras têm um carácter apenas ilustrativo, já que dispomos de resultados analíticos. 
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coordenada espacial 
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Figura 5.2: Derivada espacial do campo de velocidade para condição inicial 
periódica. 



lim u {x,t) = u±, (5.42) 
onde u± é o valor da descontinuidade. Pela conservação do momento total, vem: 



d f , , , f^i') , du , du dXit) ^ , , , 

O = — / dxu (x, t) = / dx— + / dx— H (m_ - u+) <^ 5.43) 

dtJ ^ ^ J-oo dt Jx{t) dt dt ^ ^' ^ ' 

1 Í^^^K Q 1, N 1 /"^ , ^ 2. X dXit) , , ,^ 

4^ - / dx—u (x, t) + - dx—u (x, t) = — («_ - «+) -^^ (5.44) 

2 j-oo OT 2 Jx(t) ox dt 

(í) 1 , , 

^—^ = -{u.+u+), (5.45) 

dando-nos a equação para o movimento do choque(fig |5.3|) . Levando a cabo um 
procedimento similar para a conservação de energia, obtemos 

f = -l (6.46) 

Mas, como no caso de viscosidade positiva tínhamos dissipação, por consistência, 
devemos ter no caso de viscosidade nula, > Esta dissipação é interpretada 
pela absorção do choque X{t), das partículas inicialmente dentro do intervalo la- 
grangeano [a_, a+]. Assim, uma vez absorvidas, estas partículas têm a mesma veloci- 
dade do choque. Agora, pela solução implícita, devemos ter uo (a±) = u±. Portanto, 
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a posição euleriana de uma partícula, inicialmente em a está definida, sendo que 
chamaremos, doravante, a transformação de coordenadas euler-lagrangeanas correta 
de Xe {a, t) que, como vimos, c igual a X (a, t), fora dos intervalos de choque, e igual 
a X (í) , no interior destes intervalos. 
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Figura 5.3: Decaimento da Solução 



5.3 A Equação Estatística 

Nesta secção, mais que nas anteriores, seguiremos os passos do artigo de Fournier 
e Frisch, já que este artigo é bastante intuitivo sem, contudo, lançar mão de idéias 
fenomenológicas espúrias. 

Vamos supor um campo inicial homogéneo e ergódico, assim {uq (a)) = O, e a 
covariância será 

U {x - x,t) ^ {u{x,t)u{x,t)) (5.47) 
e o espectro de energia que, é sua transformada de Fourier, será 

E{k,t)^^ í dxe-^^^^U (x, t) . (5.48) 

Consideremos a transformada de Fourier do campo de velocidade 

u{k^t) = ^j dxe-''''-'u {x, t) . (5.49) 
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Logo, teremos 

{k, t) u{k,t)) = ôik + k)E {k, t) . (5.50) 

Podemos introduzir as chamadas RFL (Representações Fourier-Lagrangeanas) se es- 
crevemos a transformada de Fourier do campo de velocidades, tomando x — Xg (a, t). 
Ou seja, 

ú(k,t) = ^ í dae-'^^^uo (a) ^X^. (5.51) 
Agora, por integração por partes, vêm as três RFL 

{ú{k,t)u{k,t)) = 1^0(^0^x^0(0) lio (ã)£xe^lee-<'''^+'^^)^ 

E (k, t) = ^{If J ^«e"''"' (^0 (a) úo (a) e-^^^KÍ»'*)-^^^'*)]) (5.52) 
E{k,t) = 1-1^1-^' Jdae-'''^[C{t,k,a)-C{t,k,Q)], 

onde 

^{a,t)^Xg{a,t)-a (5.53) 
é o deslocamento lagrangeano, e 

C{t,k,a)^ ^e-ífeKM)-í(o,í)] ^ _ (5 54) 

5.3.1 Condição Inicial Com Tempo de Pré-Choque Estatístico Positivo 

Muito embora cada realização tenha um tempo de pré-choque, a existência de um 
limite inferior, num sentido estatístico, para o tempo de pré-choque no caso aleatório, 
depende da medida de probabilidade sobre o espaço de probabilidade considerado. 
Em outras palavras, depende da distribuição de probabilidade do gradiente da ve- 
locidade. Por exemplo, se o campo inicial for gaussiano, há uma probabilidade 
finita, mesmo para argumentos tendendo a menos infinito, assim sendo o tempo de 
pré-choque estatístico é = 0. 

Por outro lado, podemos tomar condições iniciais tais que P [ú (a) < —M] — O, 
ou seja, o conjunto das realizações que apresenta t^, menores que l/M é nulo, de 
tal sorte que = l/M. Assim, podemos dizer que a distribuição é inferiormente 
limitada. Então, caracterizemos este caso. 
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Na fase regular (t < T^,), podemos calcular a enstropia pois, 

(t) = 1 [da = , (5.55) 

^ ^ 2 i l + tú{a) \ 1 + te (a) / ' ^ ^ 

onde utilizamos a hipótese de ergodicidade, denotando, para facilitar, (a) = w, já 
que a média não pode depender do ponto calculado, reescrevemos a igualdade acima 
como: 



oo /).>2 



W 

n(t)= I dwPiw) , (5.56) 

J-i 1 + tw 

sendo que aqui estamos assumindo que, — M = T* = 1. 

Também assumimos que a medida de probabilidade do gradiente é absoluta- 
mente contínua em relação à Medida de Lebesgue, ou seja, o conjunto dos pontos 
de choque eulerianos tem medida de probabilidade nula. Hipótese que é válida, 
no nosso entendimento, somente fase regular e para tempos nao muito 

próximos de T^,, pois se uma parte substancial das realizações tiver seu tempo de 
pré-choque próximo a T^,[|, como no caso gaussiano|]49[|, por exemplo(mais especi- 
ficamente como o movimento browniano), o conjunto dos pontos regulares (fora do 
intervalo de choque) tem dimensão Hausdorff 1/2, fenómeno este vinculado ao "au- 
mento" do número de choques. Assim, que se há um número grande de choques 
maduros num dado instante, obviamente a distribuição não pode ser absolutamente 
contínua à de Lebesgue, pois os choques trazem em si uma descontinuidade nos 
campos de velocidade e de seu gradiente. A única exceção seria se os choques se 
tornassem densos na reta e, portanto, restabelecêssemos a condição da medida de 
probabilidade ser absolutamente contínua com relação à Medida de Lebesgue, sendo 
o único detalhe que, em vez de termos a distribuição de probabilidade regular(fora 
do intervalo de choque), teríamos a probabilidade dos intervalos de choque. 

Pela segunda RFL podemos escrever a enstropia como: 

/oo 
^dkPE{k,t). (5.57) 

Portanto, o comportamento crítico de Í2 é determinado pela forma da distribuição 
P (w), na vizinhança do limite inferior da integral. Vamos supor que 

P (w) oc (w + 1)" . (5.58) 

<r*+e ^ a + O (e), onde O < e < a < 1. 



TQu seja, se P 



L ú{a) 
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A condição de normalização nos diz que a > 1 e, observando o comportamento, 
ainda na fase regular, teremos 

cte, a > O 

-ln(T, -í), a = , íTT; (5.59) 
(T, - í)", -1< a < 0. 

Portanto, a regula a acumulação dos tempos de pré-choques das realizações, logo a 
não divergência para a > O significa, em um senso amplo, que a maior parte das real- 
izações tem um tempo de pré-choque bem acima do tempo de pré-choque estatístico. 
Não obstante, a enstropia é infinita para tempos posteriores a T*. Seguindo um mo- 
delo simples, Fournier e Priscli(Apêndice B) obtiveram o seguinte comportamento 
para o espectro de energia: 

E{k,t) ^ A;-(^+2"/=^)F« [A;/iíi (í)] , A; ^ oo (5.60) 
Ki{t) ^ (T, -í)-5, íjT, (5.61) 

onde Fa é função positiva de decrescimento assintótico rápido, ou seja, nos instantes 
imediatamente anteriores ao tempo de pré-choque estatístico, o espectro assintótico 
ganha uma lei de potência, sendo que seu domínio aumenta quando t tende ao tempo 
de pré-choque estatístico. Isto nos leva a ver a como um tempo adimensional que 
regula a fração de realizações, que tem seu tempo de pré-choque ao redor do tempo 
de pré-choque estatístico. Sendo a = —1/2 a condição de todas as realizações terem 
o mesmo tempo de de pré-choque, sendo que esta fração vai se diluindo à medida 
que a aumenta, podemos observar, também, que esta lei de potência é mais forte 
do que a da fase dissipativa k~^, como veremos a seguir. 
A terceira RFL na fase dissipativa se escreve 

E [k, t) = í da ^e-*M^e(a,í)-Xe(o,í)]^ _ ^5_g2) 

Lembrando que Xe é constante dentro dos intervalos de choque, podemos conside- 
rar somente a contribuição dos eventos, de forma que a origem e h estejam dentro 
do mesmo intervalo de choque pois, caso contrário, a exponencial oscilará muito 
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rapidamente, cancelando a contribuição, já que desejamos somente informações 
assintóticas. Portanto, podemos escrever em uma primeira aproximação 

E {k, t)^I (t) A;-^ A; ^ oo, (5.63) 

onde / (í) é essencialmente o comprimento médio dos intervalos de choque, que 
depende das condições iniciais. 

5.3.2 Distribuição de probabilidade do gradiente da velocidade 

A distribuição de probabilidade para o gradiente da velocidade pode ser escrito 
como: 

P {x, t) = ('^ - ^« ^)) ) • (5-64) 

Se tomamos o fluxo como homogéneo e ergódico, podemos reescrever a equação 
acima como: 

f «■ *> = £So è £ (« - ^" *>) • ^"-'"^ 

Também, podemos escrever esta expressão na versão lagrangeana, que é 



1 fL 

P (C, t) = lim — / da 

L— >oo ZL J-L 



dXe (a, t) 



da 



(5.66) 



Mas como no intervalo de choque — O, podemos nos concentrar sobre o 

conjunto dos pontos regulares em um dado momento dentro do intervalo [—L,L], 
denotado por Rl (t), para o qual o envelope convexo é igual a função, logo podemos 
reescrever a expressão acima como: 

P(e,í)=lim^/ da\l + tú{a)\sU--^^], (5.67) 
L-^oo 2L Jrlíí) \ 1 + tu{a) J 

denotando por bk (^ = bk os números, de forma que: 

e - = 0. (5.68) 

Portanto, podemos escrever p{^,t) como: 

1 1 1 f 

P iC, t) = q 1 — 71-TT lini —r / daô (a — bk) ■ (5.69) 

' \l-t^f^\Ú{bk)\L^oo2LjRUt) ^ ^ ' 
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Se desejamos estudar o comportamento assintótico ^ — oo, é importante perceber 
que os pontos bk estarão nas proximidades dos pré-choques a^. Agora, fazendo uma 
expansão de ( p.68[ ) em torno do pré-choque mais próximo, vem: 



*\2 



ih - al) 



1 t 



(5.70) 



onde t* = —l/úo{al). Logo, ou temos raízes duplas ou nenhuma raiz, denotadas 
por b^. Considerando tempos anteriores ao tempo do k-ésimo pré-choque (^í < í*^, 
e a existência de raízes l + — j* < O, podendo estas duas desigualdades serem 
unidas na forma 

-]<úoial)<^-j^, (5.71) 
analogamente para t > t*, vem: 



t 2t^Í ^ t 



condições estas que podem ser resumidas em: 



-1 



1 



-1 



777 < Uq («1,1 < —. 



(5.72) 



(5.73) 



Seja agora II [a, í,^] a função indicador dos conjuntos de pontos que respeitam a 
desigualdade acima, dentro de Rl (t). Expandindo, u'q^ (al) em torno do pré-choque, 
podemos escrever ( |5.69|) como: 

h [a, t, ^] 



p (^, t) ^ \ \, lim 



-— / day^ 

2L jRr ít) V 



2L jRUt) ^ 1^(3) tu, {a 



(a - al) , 



que pode ser reescrito como: 

(2t)^^'^ 1 




Mn^^ (a) 



1 + ^ + íito (a) 



1/2 



(5.74) 



(5.75) 



Agora escrevemos esta expressão em termos da distribuição conjunta das três primeiras 
derivadas, tal que estejamos em um conjunto de pontos regulares sobre a reta 
P3 (ii, ií, u^3), Roo (í))- Logo vem: 

1/2 



1/2 



du^S) 



-1 _ 1 

t 7^ 



dú^ (ú,0,u^^\Roo(t)) 
- — V- P 



l+-+tú 



(5.76) 
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fazendo -u ~ vem: 

j(í)ier'/', e^-oo, (5.77) 

sendo esta uma característica do fluxo válida para todos os tempos e independente 
das condições iniciais, ou seja, uma característica da equação de Burgers. 

5.4 Observações 

Como este capítulo nos mostrou, a existência de singularidades nas soluções é o 
fator preponderante na deformação da estatística inicial, no caso da equação de 
Burgers. O mesmo tipo de lei foi provado recentemente]^], para o caso com uma 
força impulsiva no tempo e suave no espaço, mostrando que a lei dos 7/2 é um traço 
dos pré-clioques. 

Para o caso com força do tipo ruído braço no tempo e suave no espaço, a con- 
trovérsia ainda perdura, embora E. et a//[|TTl, [12| , |13| sustentem a existência de 



potência em 7/2 e Kraiclinan||21||, usando uma equação tipo Fokker-Plank, argu- 
mente em prol da potência em 3. Contudo, recentemente, o próprio KraichnanP? 
favorece a lei dos 7/2, mas sem obter qualquer resultado definitivo. 

A dificuldade está em saber se os pré-clioques estão presentes e separados, o que 
levaria à Lei dos 7/2, ou se existe alguma forma de agrupamento entre eles, o que 
levaria à uma potência em 3. 
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Conclusões 



Este trabalho teve como objetivo descrever uma parte da pesquisa atual feita em 
turbulência. Obviamente, não foi possível entrar nos detalhes de todas as técnicas 
utilizadas mas acreditamos ter sido cumprida a proposta básica do projeto. Dada a 
generalidade da tese, podemos agora fazer incursões por diversos problemas na área. 
Uma interessante linha de pesquisa é a possibilidade de utilização de redes neurais 
para solução de EDP's, já que isto poderá propiciar um método computacionalmente 
barato de fazer estatística no espaço de ía.se{ensemble) , ao invés de lançarmos mão 
da hipótese ergódica. 
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